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1.  INTRODUCTION 


The  location  of  a  remote  HP  transmitter,  starting  with  tha  measured  wava 
frequency  and  tha  azimuth  and  alavation  angle#  of  an  incoming  wava  front  at  a 
singla  site,  ia  typically  baaad  on  an  oversimplified  eharactarisation  of  tha 
ionoapharic  propagation  medium.  Thia  report  ia  concerned  with  improvement  of 
aingla  site  location  (SSL)  accuracy  through  improved  ionoapharic 
characterization  by  means  of  sounders.  Which  probe  tha  ionosphere  along  a 
relatively  small  sat  of  paths,  either  at  vertical  incidence  or  at  oblique 
incidence.  Tha  challenge  is  to  convert  these  measurements  to  a  complete 
characterization  of  the  ionosphere  along  all  raypaths  of  interest. 

At  present,  SSL  techniques  are  typieally  based  on  Figure  1,  Which  results 
from  the  Breit-Tuve  and  Martyn  equivalent  Path  theorems  in  radio  wave 
propagation  theory  [e.g.,  Budden,  1965  or  Davies,  1969].  The  raypath  in 
Figure  1  is  shown  by  the  arrowed  line  from  the  transmitter  to  the  receiver. 
The  measurement  determines  the  wave  frequency  f,  the  elevation  angle  B0  , 
and  the  azimuth  angle.  Which  determines  the  plane  of  propagation.  Sach  side 
of  the  constructed  isosceles  triangle  in  Figure  1  is  half  the  group  path 
length  P'  for  ray  propagation  at  frequency  f,  and  h*  is  the  virtual  height 
measured  by  a  vertical-incidence  sounder  at  the  equivalent  vertical  frequency 


(f  cos  i)/fc 


where  i  is  half  the  apex  angle,  as  shown  in  Figure  1.  The  validity  of  this 
figure  depends  on  the  approximation  that  the  ionosphere  is  horizontally 
stratified  for  the  portion  of  the  raypath  in  the  ionosphere,  i.e.,  the 
ionosphere  is  flat  and  without  tilts  for  this  raypath  portion.  Magnetic  field 
effects  are  also  ignored.  The  factor  k  in  (1),  Which  was  introduced  by 
[Smith,  1939]  as  a  rough  correction  for  ionospheric  curvature,  varies  from  1 
to  about  1.18  as  range  D  increases  from  0  to  4000  km.  The  central  angle  X  and 
distance  A  are  very  simply  related  to  D  by  the  geometry  of  Figure  1: 


X  -  D/(2Rq) 


A  -  Rq  (1-cos  X ) 


Determination  of  the  HF  transmitter  location  typically  proceeds  from 
Figure  1  and  measurements  on  the  incoming  wave  front  as  follows.  Virtual 
height  h*  versus  frequency  is  determined  at  some  point,  e.g.,  at  the  receiver, 
by  a  vertical-incidence  sounder.  Actually,  an  oblique-incidence  sounder  can 
also  be  used  in  conjunction  with  Figure  1  to  obtain  this  information.  It  is 
assumed  that  this  virtual  height  profile  is  the  same  throughout  the 
ionospheric  portion  of  the  raypath  in  question.  The  range  D  stay  then  be 
obtained  by  means  of  an  alteration  of  the  customary  [e.g.,  Davies,  1969] 
transmission  curve  technique,  based  on  Figure  1.  The  transmitter  location  is 
thus  achieved,  since  the  receiver  location  and  orientation  of  the  plane  of 
propagation  are  known.  Alternatively,  the  SSL  process  can  be  completely  and 
simply  computer-automated  [Reilly,  1983al.  First,  a  digital  file  consisting 
of  a  finite  number  of  points  (e.g.,  >  20)  along  the  virtual  height  profile  is 
stored  in  computer  memory.  Some  arbitrary  guess  is  then  made  for  the  range  D, 
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and  this  initiatss  computer  iteration  with  the  following  algorithm,  which  is 
baaed  on  Figure  1  and  Equations  (l)-(3) : 

(a)  Calculate  fv  from  (1)  and  h* (fv)  from  linear  interpolation 
of  the  h*  data  file. 


(b)  Determine  P'/2  from 

P*/2  -  (h*(fv)+  A) /cos  i 


Obtain  a  new  range  from 


.!  T  <p,/2)cos^ 
[rq  +  <P  *  /  2) 


X  •  tan 


sin# 


(d)  Test  for  convergence  on  D.  If  the  convergence  criterion  is 
satisfied,  output  the  transmitter  location.  If  not,  go  back 
to  (a)  and  repeat. 


It  was  found  IBeilly,  1983a]  that  convergence  was  rapid  and  dependable,  even 
in  cases  Where  the  initial  guess  for  D  was  grossly  inaccurate. 

By  far,  the  most  serious  approximation  of  the  preceding  paragraph  is  the 
way  that  ionospheric  assessment  is  handled.  The  ionosphere  changes  along  the 
raypath,  and  the  largest  improvement  in  SSL  will  come  from  taking  this  into 
account  properly.  An  opportunity  for  analysis  of  this  problem  comes  from  the 
SSL-BCT  (Single  Site  Location  Baseline  Certification  Test)  sounder  data  in 
December,  1982.  A  network  of  oblique  sounders  was  deployed,  principally  for 
purposes  of  frequency  management;  it  is  shown  in  Figure  2.  Six  oblique 
sounder  paths  (or  "zones'*)  are  shown,  associated  with  the  receiver  K  at  Fort 
Ord  and  the  six  transmitter  positions  T1-T6.  An  oblique  ionogram  can  be 
converted  into  an  ionospheric  true  height  profile.  Of  course,  the  oblique 
sounder  probes  the  ionosphere  all  along  the  sounder  path,  but  it  trill  be 
confirmed  that  the  true  height  profile  obtained  characterizes  the  ionosphere 
near  midpath.  The  midpath  points  X1-X6  for  each  sounder  path  are  shown  in 


Figure  2.  It  will  also  be  shown  how  each  true  height  profile  solution  can  be 
eccurately  characterized  by  only  a  few  parameters,  nine  in  most  cases.  This 
enables  economical  bookkeeping  of  ionospheric  density  profile  solutions 
throughout  the  day. 

An  oblique  ionogram  is  obtained  for  each  sounder  path  at  fifteen  minute 
intervals.  This  report  processes  all  available  oblique  ionograms  during  an 
arbitrarily  chosen  time  period  which  spans  1600  (UT)  on  December  6  to  0400  on 
December  7  (days  340  and  341).  lo  oblique  ionograms  are  available  between 
0400  and  1600  on  either  day.  The  period  chosen  includes  most  of  the  daylight 
hours  (local  time  is  eight  hours  earlier)  and  a  few  nighttime  hours,  during 
which  ionospheric  variation  at  points  X1-X6  is  thus  obtained.  Also  available 
for  comparison  during  this  period  is  the  data  from  four  vertical  incidence 
sounders,  one  at  K  in  Figure  2  and  three  others  at  the  positions  given  by  the 
solid  triangles  in  Figure  2.  Details  of  the  sounder  positions  are  given  in 
(Daehler,  1983].  It  will  be  noticed,  for  example,  that  the  DBS  (Desert 
Besearch  Station)  vertical- incidence  sounder  position  is  very  close  to  the 
midpath  position  for  Path  4.  and  this  is  used  to  demonstrate  that  the  true 
height  profile  solutions  obtained  from  the  oblique  ionograms  for  Path  4 
correspond  to  the  ionosphere  near  the  midpath,  even  when  ionospheric  tilts  are 
observed. 


fcvv'lsi 


The  ionospheric  density  variation  at  each  of  tha  control  points  X1-X6  will 
ba  asm  to  bo  noisy,  mainly  duo  to  tha  prasanca  of  travailing  ionospheric 
disturbances,  but  the  background  variation  with  tins  will  ba  seen  to  ba 
reasonably  smooth  and  predictable.  Indeed,  tha  variation  of  tha  ionospheric 
parameters,  whan  plotted  against  tha  variation  of  tha  cosine  of  the  solar 
zenith  angle,  will  be  seen  to  show  a  noisy  variation  about  well-defined  trend 
lines,  which  can  be  used  to  predict  background  ionospheric  properties  hours 
ahead  at  each  of  the  points  X1-X6.  This  will  be  useful  for  SSL  if  the 
information  at  these  points  can  be  interpolated  to  provide  ionospheric 
information  along  unsounded  paths.  It  trill  be  seen  that  this  is  possible 
through  straightforward  manipulation  of  the  trend  line  coefficients  for  each 
ionospheric  parameter.  The  role  of  other  ionospheric  models  will  also  be 
discussed  in  this  context.  Once  the  methodology  for  obtaining  ionospheric 
properties  along  unsounded  paths  is  developed,  modifications  for  improved  SSL 
techniques  follow,  and  these  are  discussed.  The  results  are  discussed  and 
future  directions  are  indicated  in  the  final  section. 

2.  FROM  OBLIQUE  IONOGRAMS  TO  TRUE  HEIGHT  PROFILES 

2A.  The  Basic  Method 

The  method  developed  [Reilly,  1984a]  to  convert  an  oblique  ionogram, 
which  shows  pulse  time  of  transit  between  transmitter  and  receiver  vs.  wave 
frequency,  to  a  true  height  profile,  which  is  height  vs.  plasma  frequency 
squared  (proportional  to  electron  density),  is  partially  based  on  a 
generalization  of  the  Breit-Tuve  theorem.  It  compares  favorably  with  all 
other  known  methods  for  obtaining  true  height  profiles  from  oblique 
ionograms.  It  employs  the  usual  simplifications.  Which  are  common  to  all 
these  methods: 

(a)  absolute  values  of  group  path  delay  times  are  known,  or  can  be  found, 

(b)  the  ionosphere  is  spherically  symmetric  (no  tilts),  and  (c)  magnetic  field 
effects  are  negligible. 

Figure  3  shows  a  typical  oblique  ionogram  and  a  sketch  of  a  few  raypaths 
by  Which  a  pulse  can  propagate  from  transmitter  to  receiver.  One-hop  (i.e., 
involving  one  ionospheric  reflection)  low-  and  high-ray  modes  are  shown,  as 
well  as  a  two-hop  mode.  The  ionogram  is  a  plot  of  relative  time  delay  vs. 
frequency,  since  time  synchronization  between  receiver  and  transmitter  is 
typically,  as  here,  not  sufficiently  accurate  to  obtain  absolute  time  delays. 
The  ionogram  in  Figure  3  shows  a  weak  one-hop  B  layer  mode  in  the  portion  of 
the  trace  just  below  2  msec  time  delay,  a  one-hop  F  mode  just  above  2  msec, 
and  a  two-hop  F  mode  Which  starts  near  3  msec.  In  order  to  convert  relative 
time  delays  to  absolute  time  delays,  an  assumption  is  made  about  the  height  of 
the  bottom  of  the  E  layer,  which  is  associated  with  the  zero-frequency  time 
delay  of  the  one-hop  E  mode.  In  most  cases  this  height  is  assumed  to  be  at  90 
km,  although  in  the  nighttime  case,  with  sporadic  E  evident,  this  value  is 
frequently  changed  to  100  km.  It  will  be  seen  that  simplification  (a)  above 
can  thus  be  sufficiently  well  justified. 

The  simplification  (b)  above  can  be  regarded  as  the  first  step  of  an 
iterative  process  in  which  ionospheric  density  profiles  are  inferred  at  the 
various  sounder  control  points,  ionospheric  tilts  are  inferred  from  this 
information,  and  then  ionograms  are  reduced  from  this  tilt  information  in  a 
hypothetical  second  iteration.  It  will,  however,  be  argued  that  the  ionograms 
really  represent  the  effect  of  the  ionospheres  near  the  midpath  points  in  the 
presence  of  actual  tilts,  and  that  simplification  (b)  is  therefore 
sufficiently  justified  for  our  purposes.  In  other  words,  the  ionogram 
wouldn't  change  very  much  if  the  midpath  ionosphere  replaced  the  actual  one 
throughout  the  sounder  path. 
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Simplif ication  (e)  is  justified  for  obliqus  ionograns,  where  magnetic 
fiald  ef facta  only  amount  to  a  relatively  snail  splitting  of  tha  traca  in  tha 
vicinity  of  tha  MUF  and  up  into  tha  high-ray  portion  of  tha  traca.  This  is 
saan  in  Figure  3.  It  is  satisfactory  to  astinata  tha  no-fiald  traca  valuas  by 
choosing  points  about  halfway  between  tha  split  ordinary  node  and 
extraordinary  node  traces  [Kopka  and  Holler,  1968]. 

Tha  procedure  adopted  to  process  ionograns,  such  as  the  one  in  Figure  3 
fron  the  SSL-BCT  [Daehler,  1983],  is  to  first  identify  the  one-hop  nodes,  and 
then  enter  coordinates  sequentially  fron  the  one-hop  trace  into  a  BASIC 
language  conputer  progran,  written  for  tha  Tektronix  4052  microconputer.  The 
interface  between  tha  ionogran  and  tha  progran  is  a  Tektronix  4956  Graphics 
Tablet  Digitiser.  The  one-hop  trace  starts  at  the  bottom  of  tha  B  layer 
(estinated  at  f «0)  and  advances  up  in  reflection  height  to  the  topmost 
observed  portion  of  the  F  layer.  As  in  Figure  3,  tha  trace  is  typically 
discontinuous,  showing  a  portion  of  the  B  layer  low-ray  trace  (e.g.,  below  the 
B  MUF),  a  portion  of  the  F  laye?  low-ray  trace,  and  ending  after  a  relatively 
small  portion  of  tha  F  layer  high-ray  trace  (e.g.,  above  the  F  MUF).  As 
points  are  entered  fron  the  one-hop  trace,  tha  progran  automatically  obtains 
true  height  profile  solutions.  Tha  ionospheric  density  profile  is  thus  built 
up  fron  tha  bottom.  Details  of  the  method  are  contained  in  (Beilly,  1984a] . 
Examples  will  be  shown  below. 
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2B. _ Idealized  True  Height  Profile. Fit 

It  is  convenient  to  define  a  model  density  profile,  aa  shown  in 
Figure  4.  It  is  composed  of  a  series  of  segments ,  in  each  of  Which 
R2fg  2  varies  either  linearly  or  quadratically  with  R,  Where  R  is  tha 
radial  distance  fron  the  center  of  tha  earth.  Tha  height  h  is  found  by 
subtracting  the  radius  of  tha  earth  fron  R.  A  dependence  of  this  forn,  i.e.. 
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where  A,  B,  and  C  are  constants,  enables  tha  integrals  for  group  path  delay 
and  range  to  be  evaluated  in  closed  forn  [Reilly,  1983a],  and  this  has  been 
found  convenient  for,  e.g. ,  the  computation  of  vertical  and  oblique 
ionograns.  The  idealized  profile  in  Figure  4  consists  of  a  linear  segment  for 
the  B  region,  a  linear  segment  for  the  B-F  transition  region,  a  linear  segment 
for  the  lower  part  of  the  F  region,  and  a  parabolic  segment  for  the  upper  part 
of  the  F  region  up  to  the  F  layer  naxiraun.  Bote  that  the  segments  of 
R2fg  2  are  linear  or  parabolic,  but  the  segments  of  fg  2  also  appear 
linear  or  parabolic,  as  in  Figure  4,  because  R  is  nearly  constant  in  the 
altitude  range  shown.  The  profile  is  specified  by  nine  parameters,  the  five 
height  parameters  hbE,  hmB,  hbF,  hmFl,  and  hmF2 ,  and  four  density  parameters 
fOB2,  fbF2,  fOPl2,  and  fOF22.  The  semi-thickness  parameter  Y2  for  the 
parabolic  layer  is  redundant.  These  are  suggestive  labels  only,  not 
necessarily  having  their  usual  ionospheric  significance.  For  example,  the  B 
layer  maximum  is  not  to  be  found  at  (fOE2,  hmB),  nor  does  (fOFl2,  hmFl) 
necessarily  give  the  coordinates  of  the  FI  layer  maximum.  Instead,  these 
parameters  are  manipulated  by  the  computer  progran  to  give  the  best 
least-squares  fit  of  the  idealized  profile  in  Figure  4  to  the  computed  true 
height  profile.  The  latter  is  defined  by  the  set  of  calculated  solution 
points,  which  are  interpolated  by  segments  in  each  of  which  the  dependence  of 
R2fg  2  on  R2  is  linear  [Reilly,  1984a].  Indeed,  the  fit  of  the 
idealized  profile  is  so  good  that  the  calculated  true  height  profile  is  well 
specified  by  these  nine  parameters. 
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The  variation  of  g  =*  R2fg  2  in  each  of  the  linear  segments  of  the 
idealized  profile,  in  Which  g  varies  from  gj_  to  g2  as  h  varies  from  h^ 
to  h2,  is 

8  “8 

g  -  gj  +  ~~r  h|  ■  (h-hj)  (linear)  (5) 

The  variation  in  the  parabolic  segment,  in  Which  g  varies  from  gOFl  to  gOF2  as 
h  varies  from  hmFl  to  hmF2,  is 

g  -  gOF2  £l-  ^  hmF2-hmFl  )  (  1-  gOFl)  j  (parabolic)  (6) 

An  example  of  the  computer  solution  process  is  shown  in  Figures  5-7. 
One-hop  trace  points  are  fed  into  the  computer  program  from  the  ionogram  of 
Figure  3,  and  the  resulting  true  height  profile  solution  points  are  shown  by 
the  pluses  in  Figure  5,  connected  by  solid  lines.  The  program  would  generate 
a  monotonically  increasing  profile  (fg  2  increasing  uniformly  with 
altitude),  except  for  the  presence  of  a  reentrant  procedure  [Reilly,  1984a], 
which  may  be  summarized  as  follows.  Whenever  (except  in  the  E  layer)  three 
adjacent  monotonic  solution  points,  n-1,  n,  and  n+1 ,  are  such  that  the  segment 
between  n  and  n+1  solution  points  has  a  slope  d(R2fg  2)/dR2  twenty 
percent  greater  than  that  for  the  segment  between  points  n-1  and  n,  the 
reentrant  procedure  is  initiated.  A  new  reentrant  point  is  inserted  after 
point  n-1,  placed  on  the  line  upward  from  n-1  with  slope  dh/dfg  2  * 

-  /z. s  km/MH z2  at  an  altitude  such  that  the  segment  from  this  reentrant  point 
to  the  recomputed  solution  point  n  has  the  same  slope  as  the  segment  from  n  to 
the  recomputed  solution  point  n+1.  If  and  When  this  condition  is  satisfied, 
the  new  point  and  the  latter  two  points  are  redefined  as  points  n,  n+1,  and 
n+2  in  the  solution  set.  The  solution  procedure  then  continues  on  to  the  next 
point  of  the  one-hop  trace.  Evidently,  the  reentrant  procedure  was  carried 
out  several  times  in  Figure  5  between  the  E  and  F  layers. 

After  the  last  point  on  the  one-hop  trace  has  been  converted  to  a  true 
height  profile  solution  point,  the  program  finds  the  best  parabolic 
least-squares  fit  of  the  last  five  solution  points  and  determines  the  F2  peak 
parameters  (fOF22,  hmF2)  from  this  parabola.  The  extrapolation  of  the 
parabola  upward  and  backward  from  the  last  solution  point  is  shown  by  the 
dotted  line  in  Figure  5. 

The  program  next  accepts  user  input  for  the  designation  of  (hmE,  fOE2) 
as  one  of  the  solution  points  (J)  and  an  initial  guess  for  the  parameters  hbF 
and  hmFl  (HI  and  H2).  The  parameters  hbF,  hmFl,  fbF2,  and  fOFl2  are  then 
determined  by  the  program  in  an  iterative  procedure  which  finds  the  best 
least-squares  fit,  consistent  with  continuity  requirements,  of  the  calculated 
true  height  profile  by  the  idealized  profile  between  the  (hmE,  fOE2) 
solution  point  and  the  fifth  to  last  solution  point.  Beyond  this  latter  point 
the  idealized  solution  is  the  aforementioned  parabolic  fit.  The  resulting 
best  idealized  fit  is  shown  by  the  dotted  line  in  Figure  6.  It  is  good  enough 
to  characterize  the  calculated  true  height  profile,  which  is  thus  equivalent 
to  the  nine  parameters  given  in  Figure  7  (Y2  is  redundant) .  This  is  a  great 
bookkeeping  simplification. 
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Calculated  true  height  profile  points  (pulses)  and  best  quadratic  least  squares  Tit 
to  last  five  points  (dotted) 
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Calculated  profile  and  best  least  squares  fit  of  idealized  profile  (dotted) 
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It  appears  that  ttea  uaar  can  almost  entirely  ba  eliminated  from  tha 
solution  proeass,  as  follows.  Tha  ana-hop  traca  can  ba  downloadad  to  tha 
computer  program  froa  tha  souadar  as  a  digital  fila.  This  would  aaka  it 
uraiaeassary  for  a  usar  to  march  up  tha  ona-hop  traca  of  tha  ionograa  with  his 
digitising  davica.  Tha  solution  point  (J)  for  (hraK,  fOE2)  in  Figure  S  could 
ba  datarminad  by  tha  eoaputar  as  tha  point  just  bafora  tha  first  rasntrant 
point  or  tha  last  point  balow  h-130  km,  whichavar  comas  first  (racall  that  tha 
reentrant  procedure  is  not  permitted  until  ona-hop  traca  points  in  tha  F  layer 
are  processed).  Tha  initial  guess  for  hbF  and  hmFl  in  Figura  5  could  be  made 
by  tha  computer  as  tha  last  reentrant  point  height  and  a  value  10  km  lass  than 
tha  fifth  to  last  solution  point  height,  respectively.  These  and  other  steps 
make  tha  prospects  bright  for  complete  computer  automation  for  tha  vast 
majority  of  ionograms  encountered. 

A  condition  can  occur,  especially  at  nighttime,  whan  eleven  parameters  are  ' 
suitable  to  describe  the  true  height  profile  solution.  In  this  case,  tha 
trial  altitude  of  the  first  reentrant  point  advances  so  far  up  tha  line  with 
slope  -/2.5’lon/MHz2  that  it  reaches  tha  fg  2*0  axis  bafora  tha  reentrant 
condition  is  satisfied.  This  point  of  intersection  with  tha  axis  becomes  tha 
definition  of  (fOB2,  hmS).  From  hare,  the  trial  reentrant  point  advances  up 
the  axis  to  a  point  such  that  tha  computed  slopes  of  tha  next  two  segments  are 
equal,  which  is  the  reentrant  condition.  The  fit  of  idealized  profile  to  tha 
calculated  true  height  profile  proceeds  as  bafora.  An  example  of  this  is 
shdwn  in  the  true  height  profile  solution  of  Figure  8.  Two  new  parameters  are 
defined  to  characterize  the  true  height  profile,  i.e.,  hsK  and  fsE2.  These 
and  the  other  parameters  are  indicated  in  Figure  8.  Although  eleven  are 
shown,  there  are  really  ten  Independent  parameters  in  our  present  procedure, 
since,  for  example,  hsE  is  determined  from  fsE2,  hmK,  fOE2,  and  the  fixed 
slope  value  -17. .5  for  the  reentrant  segments.  The  value  of  this  slope  for  the 
first  reentrant  point  might  at  some  future  time  be  made  flexible.  The  best 
fit  in  Figure  8  was  obtained  by  the  computer  for  fbF2«0,  although  this  is 
not  always  the  case. 

The  procedure  for  nighttime  ionograms,  such  as  the  one  which  yielded 
Figure  8,  which  shows  a  sporadic  E  layer,  will  be  to  innut  a  point  on  the 
one-hop  E  layer  trace  a  small  distance  out  from  the  fg  2«0  axis,  and  then 
proceed  up  to  the  F  layer  part  of  the  one-hop  trace.  This  defines  a  small, 
token  E  layer,  as  in  Figure  8,  which  has  no  significant  effect  on  the 
calculated  F  layer  solutions. 

There  will  be  a  few  instances  in  Which  the  first  reentrant  point  in  Figure 
8  will  occur  close  to,  but  not  on  tha  fg  2*0  axis.  Then  the 
eleven-parameter  fit  of  the  type  of  Figure  8,  except  that  fOB2  is  not  zero, 
is  appropriate  to  describe  the  calculated  true  height  profile. 

2C.  SSL-BCT  Prototype  Calculations 

The  purpose  of  this  section  is  to  test  the  solution  procedure  for 
obtaining  true  height  profiles  on  known  cases,  Which  resemble  the  conditions 
of  the  SSL-BCT.  From  the  known  sounder  geometry  in  Figure  2  [Oaehler,  1983] 
it  is  found  that  the  transmitter-receiver  separation  ranges  for  Paths  1,2, — ,6 
are  385.3,  801.0,  918.3,  1501.9,  600.5,  and  1168.3  km.,  respectively. 

Known-case,  one-hop  oblique  ionograms  are  calculated  from  the  prototype 
daytime  profile  in  Figures  9  and  10  for  the  five  distinct  ranges  from  the 
preceding  list.  In  these  figures  the  variation  of  R2fg  2  with  radius  R 
is  parabolic  in  the  E  layer  from  h*hbE,  past  its  maximum  at  (fOB2,  hmB) ,  to 
the  point  h  where  R2fg  2«RmB2f0E2Kv  2.  The  linear  E-F 
transition  region  starts  here  and  ends  at  (fbF2,  hbF).  The  next  linear 
segment  for  the  lower  part  of  the  F  layer  proceeds  from  here  to  (fOFl2, 


taPl),  and  tha  parabolic  F2  ration  extends  from  bare  to  just  past  its  maxiaum 
at  (fOP22,  hmF2) .  Thasa  parameters  for  tha  profile  in  Figure  10  are  given 
in  Figure  9,  as  are  tha  coefficients  in  lq.  (4)  for  tha  associated  five 
segments,  starting  with  tha  ground  up  to  tha  bottom  of  tha  K- layer .  K  is  in 
units  of  km,  with  tha  earth's  radius  being  6371.2  tea.  This  profile  is  similar 
to  that  defined  in  Figure  4  for  the  idealised  profile  fit,  except  that  the  K 
layer  is  linear  there. 

Figure  11  depicts  the  vertical  ionogram  associated  with  the  given  density 
profile,  and  the  next  ten  figures  (12-21)  show  the  calculated  one-hop  oblique 
ionograas  for  five  distinct  ranges  froa  the  above  SSL-BCT  list.  These 
ionograas  can  be  regarded  as  exact,  sines  they  are  based  on  an  exact 
evaluation  of  the  associated  group  path  and  range  integrals,  made  possible  by 
the  fora  of  8q.  (4).  Each  range  value  is  associated  with  a  pair  of  figures, 
the  first  of  which  gives  the  reflection  heights  and  elevation  angles  for  a  set 
of  points  on  the  one-hop  trace.  The  second  one  shows  circled  points  Which  one 
might  at  random  attempt  to  digitize  on  an  actual  ionogram  for  input  to  the 
computer  program  which  finds  the  true  height  profile  solution.  The  only  part 
of  the  one-hop  trace  that  shows  up  in  a  typical  ionogram  is  a  portion  of  the 
low-ray  B  layer  (below  the  B  MUF),  most  of  the  low-ray  F  layer,  and  a  fairly 
small  portion  of  the  high-ray  F  layer  (above  the  F  MUF).  The  circled  points 
shown  are  located  on  these  portions  of  the  trace. 

The  true  height  profile  solution  program  has  been  tested  on  the  oblique 
ionograms  in  Figures  13,15,17,19,  and  21.  Two  methods  are  eaployed.  In  one 
method  the  calculated  values  of  absolute  time  delay  vs.  frequency  for  the 
circled  ionogram  points  are  stored  in  a  data  file.  These  one-hop  trace  values 
are  read  sequentially  by  the  program,  which  generates  the  corresponding  true 
height  profile  solution  values.  The  other  method.  Which  is  used  for  a 
measured  ionogram,  requires  an  initial  specification  of  hbB,  Which  is  used  to 
convert  the  observed  relative  time  delays  to  absolute  time  delays.  The 
ionogram  is  interfaced  with  the  computer  program  through  a  magnetic  wire  grid 
tablet  device,  on  Which  the  ionograa  is  laid.  A  cross-hairs  device  is 
advanced  up  the  one-hop  trace,  and  values  of  relative  tine  delay  vs.  frequency 
are  thus  digitized  for  input  to  the  computer  program,  Which  converts  relative 
time  delays  to  absolute  time  delays  and  generates  the  true  height  profile 
solution.  The  technique  for  carrying  out  the  second  method  has  been  varied 
and  compared  with  the  first  method,  and  the  results  of  this  experimentation 
are  given  in  Table  1.  It  is  seen  that  the  solution  process  has  been  carried 
out  independently  forty  times  on  the  five  ionograms.  Bach  solution  is 
characterized  by  nine  parameters,  as  discussed  in  Section  2B.  Of  course,  it 
is  difficult  to  visualize  the  results  from  Table  1,  and  so  each  of  the  forty 
cases  is  also  plotted  in  Appendix  K,  where  the  calculated  9-parameter  true 
height  profile  of  Table  1  is  depicted  by  the  x-marks,  and  the  exact  model 
profile  of  Figure  10  is  shown  by  the  solid  line  for  comparison. 

Bach  nine-paramater  result  in  Table  1  is  followed  by  a  short  description 
of  the  points  included  for  processing  from  the  one-hop  ionogram  trace.  In 
most  cases  the  circled  points  in  Figures  13,15,17,19,  and  21  are  the  ones 
processed.  These  involve  only  the  low-ray  portion  of  the  B  layer,  in  keeping 
with  typical  data.  Inspection  of  Appendix  K  shows  that  the  B  layer  is  thus 
clipped  near  its  peak.  Inclusion  of  high-ray  B  layer  points  improves  the  B 
layer  peak  description  somewhat  (cf.  los.  8,10,11,25,  and  33  in  Table  1),  but 
the  F  layer  solution  is  fortunately  insensitive  to  this.  In  the  tablet  mode 
of  operation,  hbB  is  an  input  parameter.  In  practice  it  is  usually  estimated 
to  be  90  km  for  daytime  ionograms  and  100  km  for  nighttime  ionograms,  when 
sporadic  B  conditions  are  evident.  The  effect  of  varying  hbB  in  Bos. 
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Fig.  13  -  SSL-BCT  prototype  one-hop  ionogram 
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Fig.  15  -  SSL-BCT  prototype  one-hop  ionogram 


Fig.  16  —  SSL-BCT  prototype  one-hop  ionogram 
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Fig.  17  -  SSL-BCT  prototype 
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Fig.  21  —  SSL-BCT  prototype  one-hop  ionogram 


calculated  from  the  SSL-BCT  prototype  one-hop  ionograms 
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6,12,13,19,20,28,  29,  and  35-40  can  ba  saan  by  comparing  thaaa  raaulta  with 
tha  hbE«90  ka  raaulta  and  with  aaeh  other.  Tha  B-layar  aolution  ahowa  about 
the  expaetad  variation  with  hbB.  Soaw  additional  variability  in  tha  B  layar 
aolution,  aapaeially  at  larga  ranges,  darivaa  from  human  arror  in  tablat 
digitization  of  tha  flat  low-ray  B  layar  part  of  tha  ionograa  traca.  Tha 
affaet  of  tha  variability  of  tha  B  layar  aolution  on  tha  F  layar  aolution  ia 
found  to  ba  ralativaly  small.  Moat  SSL  moasuramonta  ara  aaaoeiatad  with  F 
layar  raflaetiona,  so  that  this  insensitivity  of  tha  F  layar  solution  to  tha 
errors  mads  for  tha  B  layar  is  fortunate  for  HF  transmitter  location.  It  is 
also  fortunate  that  tha  true  height  profile  calculations  ara  reproducible  and 
capable  of  high  accuracy,  as  saan  in  Appendix  K,  and  that  tha  tablat 
digitization  method  ia  almost  as  accurate  as  tha  data  file  method  (compare 
Mos.  1  and  4,  11  and  16,  18  and  21,  and  26  and  30).  Tha  tablat  digitization 
tends  to  give  fOF2*  values  slightly  too  large,  a  result  of  nisei ignamnt  of 
wires  in  the  grid,  but  the  effect  is  insigificant  (and  largely  correct ibla ) . 

Tha  cases  for  range  1501.9  in  Table  1  show  a  processing  difficulty  fox  tha 
B  layar  trace  that  seems  characteristic  of  larga  ranges.  It  was  found  that 
circled  points  2-7  in  Figure  21  could  not  ba  processed  by  tha  computer 
algorithm,  for  reasons  presently  not  clear.  Tha  procedure  adopted  in  the 
tablet  operation  was  to  (1)  lay  tha  cross-hairs  on  and  digitize  tha  estimated 
f«0  position  for  the  bottom  of  tha  one-hop  traca,  (2)  input  tha  associated 
value  of  hbB  (usually  estimated  to  be  90  ka),  which  was  used  to  convert 
observed  relative  time  delays  to  absolute  time  delays,  (3)  attempt  to  digitize 
the  observed  point  at  the  B  MUF,  and,  failing  a  solution  here,  (4)  advance  tha 
cross-hairs  up  the  high-ray  portion  of  the  B  layer  one-hop  trace  until  the 
second  solution  point  is  found.  This  is  the  procedure  also  followed  for  an 
actual  ionogram  at  this  range,  except  that  the  high-ray  portion  of  tha  B  layar 
trace  has  to  be  imagined,  using  Figure  20  or  21  as  a  guide.  From  hare,  tha 
points  in  the  F  layer  part  of  tha  traca  ara  processed,  as  par  usual.  Despite 
inaccuracy  in  the  B  layer  solution,  tha  F  layar  solution  is  again  quite 
accurate,  as  seen  in  Bos.  34-39  in  Table  1  and  Appendix  K. 

Most  cases  in  Table  1  employ  tha  reentrant  (KB)  procedure  (cf.  Section 
2B) ,  but  some  do  not.  There  is  little  difference  in  this  case,  since  the 
starting  profile  in  Figures  9  and  10  is  not  essentially  reentrant,  i.e.,  it  is 
almost  a  monotonically  increasing  function.  For  profiles  Which  have  a 
significant  reentrant  portion,  the  KB  procedure  usually  gives  much  greater 
accuracy  (Reilly,  1984a].  Nevertheless ,  the  RE  procedure  gives  small 
improvement  in  this  case  too.  The  solutions  of  Table  1  tend  to  underestimate 
the  value  of  hbF,  partially  because  of  the  discontinuous  slope  at  tha  boundary 
of  tha  B-F  transition  segment  and  tha  F  layer  segment  in  Figure  10.  Tha  RE 
procedure  tends  to  improve  the  solution  somewhat  in  the  vicinity  of  this 
boundary  (Reilly,  1984a). 

The  preceding  boundary  point  is  associated  with  the  cusp  at  the  low 
frequency  end  of  the  F  layer  low-ray  portion  of  the  ionogram  trace.  It  is 
important ,  at  least  for  this  model  profile,  to  select  a  point  near  the  cusp 
for  the  true  height  profile  calculation.  In  most  of  tha  range  385.3  entries 
of  Table  1,  tha  sixth  circled  point  of  Figure  13  was  selected  for  processing. 
This  is  apparently  not  near  enough  to  the  cusp.  Tha  affect  of  this  in  Mos. 

1-5  and  7-9  is  saan  in  Appendix  K  to  place  tha  calculated  solution  points  too 
low  in  altitude  for  tha  lower  part  of  tha  F  layar.  Inspection  of  Mo.  10  in 
Appendix  K  shows  tha  improvement  from  selecting  a  point  closer  to  the  cusp. 

Tha  entries  for  larger  ranges  in  Table  1  use  points  closer  to  tha  cusps  on  tha 
ionogram  traces,  and  batter  accuracy  is  thus  achieved.  This  point  is  further 
illustrated  by  comparison  with  Mos.  15,23,25,33,  and  40.  In  practice, 
however,  difficulty  in  getting  a  solution  is  sometimes  found  by  trying  to 


select  a  point  too  naar  tha  cusp.  In  those  eases,  the  digitization 
cross-hairs  are  gradually  moved  away  from  the  cusp  until  a  solution  can  he 
obtained . 

One  further  problem  with  the  computer  algorithm  should  be  noted.  If 
points  too  far  up  in  the  high-ray  portions  of  the  B  and  F  layer  portions  of 
the  ionogram  trace  are  processed,  the  associated  true  height  solutions  are 
found  to  be  too  high  {Reilly,  1984a] .  This  leads  to  an  inaccurate 
determination  of  hmF2  and,  to  a  lesser  extent,  of  fOF22.  In  practice,  the 
computer  algorithm  has  difficulty  generating  a  solution  in  some  of  these 
cases,  normally,  this  is  not  a  problem  for  typical  ionograms,  where  enough  of 
the  high- ray  trace  is  visible  to  find  the  peak  parameters  with  sufficient 
accuracy,  but  not  so  much  that  inaccurate  heights  are  calculated.  By  using  a 
little  sophistication  in  signal  processing,  it  is  often  possible  to  generate 
ionograms  that  show  ntoo  much"  of  the  high-ray  F  layer  trace.  It  is  not 
difficult,  however,  to  stop  processing  the  high-ray  points  at  some  point. 

This  is  one  of  the  values  of  prototype  calculations  of  this  type.  Which  serve 
as  a  guide  for  the  real  thing.  It  is  found  that  stopping  the  calculations  at 
or  slightly  below  the  last  circled  points  in  Figures  13,15,17,19,  and  21 
yields  sufficient  accuracy  for  the  peak  parameters  of  the  F2  layer.  In  Nos. 
10,17,25,33,  and  40  in  Table  1,  an  attempt  was  made  to  digitize  points  well 
beyond  the  last  circled  points  of  the  aforementioned  figures.  It  is  found 
that  there  are  gross  errors  made  in  hmF2  for  the  two  largest  ranges ; 
apparently  the  errors  incurred  are  not  serious  for  the  smaller  ranges.  It 
will  also  be  seen  that  the  values  for  hmFl  and  fOFl2  are  anomalously  high 
for  the  two  largest  range  cases.  This  does  not  mean  that  the  middle  part  of 
the  F  layer  solution  is  also  very  inaccurate.  Inspection  of  Nos.  33  and  40  in 
Appendix  K  shows  this.  The  significance  of  hmFl  and  fOFl2  as  parameters  in 
a  least-squares  fitting  scheme  should  be  remembered. 


2D.  Path  4  True  Height  Profiles  and  Midpath  Ionospheres 

An  oblique-incidence  sounder  senses  the  ionosphere  all  along  the 
raypaths,  which  brings  into  question  the  significance  of  the  true  height 
profiles  calculated  from  the  oblique  ionograms.  Does  such  a  true  height 
profile  represent  some  freighted  average  of  ionospheres  all  along  the  sounder 
path?  It  is  reasonable  to  suppose  that  this  is  the  case,  and  it  may  also  be 
argued  that  the  weighting  factor  for  the  ionosphere  near  midpath  is 
particularly  large,  since  this  is  where  most  of  the  refraction  of  the  rays 
occurs  for  one-hop  modes,  at  least  for  the  relatively  small  ionospheric  tilts 
(  <5°)  actually  found.  Rephrased,  the  true  height  profile  is  the  solution 
of  an  inverse  scattering  problem,  where  most  of  the  scattering  takes  place 
near  the  midpath  latitude  and  longitude  coordinates.  This  notion  is  supported 
by  analysis  of  oblique-backscatter  records  [Basler  and  Scott,  1973],  where  the 
oblique-backscatter  ionograms  are  converted  to  equivalent  vertical  ionograms, 
using  approximations  similar  to  those  discussed  for  Figure  1,  and  these  are 
found  to  correspond  quite  closely  to  vertical  ionograms  measured  at  midpath. 

The  sounder  configuration  for  Path  4  in  Figure  2  allows  a  further 
examination  of  the  correspondence  between  the  true  height  profile  solutions 
and  the  associated  midpath  ionosphere.  The  approach  is  to  identify  cases 
where  ionogram  data  from  the  vertical-incidence  sounders  at  Fort  Ord  (R  in 
Figure  2)  and  Desert  Research  Station  (DRS  in  Figure  2)  can  be  compared  with 
equivalent  data  obtained  from  the  true  height  profile  solutions  calculated 
from  Path  4  oblique  ionograms.  This  will  also  provide  another  check  on  the 
accuracy  of  the  calculated  true  height  profiles.  It  would  be  nice  to  also 
have  vertical  ionogram  data  from  Boulder,  CO,  (i.e.,  near  T4  in  Figure  2),  but 
it  seems  that  this  data  is  unreliable  during  the  time  periods  of  interest  in 


«**  V"  h*  V* 


the  SSL-BCT  data  basa.  Hanca,  during  tiaa  pariods  in  which  thara  ara  obliqua 
ionograms  for  Path  4  and  vertical  ionograms  fro*  Fort  Ord  and  Dasart  Raaaarch 
Station,  true  height  profiles  ara  calculated,  and  than  used  to  calculate 
vertical  ionograms,  as  in  Figure  11,  for  comparison  with  the  Measured  vertical 
ionograms.  The  comparison  is  properly  made  with  the  ordinary-mode  traces  of 
these  measured  vertical  ionograms,  since  the  reflection  conditions  are 
identical  between  the  ordinary-mode  and  no-field  cases.  If  the  calculated 
true  height  profiles  are  to  correspond  closely  to  the  midpath  ionosphere,  then 
the  calculated  vertical  ionogram  should  be  much  closer  to  the  Desert  Research 
Station  data  than  to  the  Fort  Ord  data  in  cases  Where  these  data  differ  from 
each  other. 

The  time  periods  chosen  for  the  above  comparison  were  December  14  (Day 
348)  between  1600  and  1915  (UT) ,  December  16  (Day  350)  between  0100  and  0245, 
and  December  18  (Day  352)  between  1530  and  1545  (two  oblique  ionograms). 
Surprisingly,  this  represents  a  significant  fraction  of  the  cases  Which  were 
available  for  this  comparison.  During  the  above  three  time  periods,  the 
ionospheres  at  DRS  and  Fort  Ord  are  expected  to  differ.  The  first  and  third 
time  periods  include  early  morning  hours  near  sunrise,  and  the  second  time 
period  includes  dusk  and  early  nighttime  hours  (subtract  eight  hours  from  UT 
to  get  local  time) . 

The  aforementioned  comparison  is  shorn  in  Appendix  A  for  seventeen 
different  times  during  the  three  time  periods.  At  each  time  a  copy  is  shown 
of  the  oblique  ionogram  from  Which  the  true  height  profile  was  calculated. 

This  is  followed  by  the  comparison  between  the  equivalent  no-field  vertical 
ionogram  (pluses  and  dots),  which  was  calculated  from  the  true  height  profile 
solution,  and  the  ordinary-mode  traces  from  the  corresponding  vertical 
ionograms  at  Desert  Research  Station  (x-marks)  and  Fort  Ord  (little  black 
rectangles) .  In  some  cases  only  one  of  the  two  oblique  ionograms  shown  on  a 
page  is  processed;  the  other  one  is  not,  because  of  an  incomplete  set  of  data 
for  comparison  at  that  time.  In  some  cases  a  line  will  be  seen  to  be  drawn 
through  the  E  layer  portion  of  the  one-hop  trace.  This  was  used  to  help 
estimate  the  position  of  the  bottom  of  the  trace  at  f=0.  in  many  cases  this  E 
layer  portion  was  weak  and  indistinct,  thus  making  it  even  more  difficult  to 
process  the  B  layer  portion  of  the  one-hop  trace  for  Path  4,  as  discussed  in 
Section  2C.  In  some  cases,  it  was  difficult  to  identify  the  ordinary  mode  on 
the  Fort  Ord  and  DRS  vertical  ionograms,  particularly  in  the  jumble  of  points 
near  the  E  layer  critical  frequency.  Several  points  of  the  jumble  were 
typically  plotted,  as  were  some  sporadic  E-layer  points  (e.g.,  see  Day  352, 
1545). 

It  is  generally  seen  in  the  vertical  ionogram  comparisons  that,  in  cases 
Where  the  DRS  and  Fort  Ord  data  differ,  the  calculated  equivalent  vertical 
ionogram  points  for  Path  4  agree  much  more  closely  with  the  DRS  data.  This  is 
in  keeping  with  the  notion  that  the  calculated  true  height  profile  solution 
characterizes  the  midpath  ionosphere.  An  interesting  exception  is  the 
comparison  shown  for  Day  350,  0200,  where  the  agreement  with  the  Fort  Ord  data 
seems  better  near  the  F2  critical  frequency.  Still,  it  must  be  noticed  that 
the  DRS  and  Fort  Ord  data  don't  differ  very  much  from  each  other.  Further, 
the  01  calculations  are  based  on  a  relatively  incomplete  ionogram.  Hence, 
this  apparent  exception  does  not  seem  at  all  significant. 

In  a  few  cases,  the  true  height  profile  solution  must  be  somewhat 
inaccurate.  Examples  of  this  are  shown  for  Day  348  at  times  1615  and  1830. 
Evidently,  the  calculated  E  layer  and  F2  layer  critical  frequencies  are  too 
high  in  each  of  these  cases.  This  is  not  surprising  for  the  E  layer  critical 
frequency,  because  of  the  highly  approximate  way  that  the  E  layer  has  to  be 
treated  for  Path  4,  as  previously  discussed.  The  fact  that  the  calculation  of 
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the  vertical  ionogram  near  the  F2  critical  frequency  differs  from  the  DRS  data 
in  some  cases  may  be  a  reflection  of  the  uncertainty  in  the  extrapolation 
procedure  used  to  infer  F2  peak  parameters  in  the  true  height  profile 
calculation.  Overall,  however,  the  agreement  between  the  DRS  data  and  the 
calculated  equivalent  vertical  ionogram  is  striking,  especially  in  the  F 
layer,  and  this  lends  confidence  in  the  calculations.  Greater  accuracy  should 
be  obtained  by  repeating  the  calculations  on  ionograms  which  are  replotted  for 
lower  threshold  values.  All  signals  are  excluded  from  the  ionogram  plots  for 
intensities  below  this  threshold,  and,  if  it  is  lowered,  greater  detail  in  the 
B  layer  and  on  the  high-ray  F  layer  portion  of  the  ionogram  trace  will  be 
visible. 

3.  SSL-BCT  CALCULATIONS  FOR  DECEMBER  6-7: 

ia  Ionograms 

The  method  developed  for  obtaining  true  height  profiles  from  oblique 
ionograms  has  been  applied  to  the  ionograms  of  SSL-BCT  for  December  6-7  (Days 
340  and  341).  From  the  discussion  of  Section  2D,  these  profiles  will  specify 
the  ionosphere  at  the  midpath  points  X1-X6  in  Figure  2  during  this  time 
period.  Later,  these  results  will  be  compared  with  available  vertical 
ionogram  data  (Daehler,  1983]  associated  with  the  stations  denoted  by  the 
solid  triangles  and  R  in  Figure  2.  The  locations  of  all  these  points  are 
given  in  Table  2. 

The  time  period  covered  is  between  1600  UT  on  December  6  (Day  340)  and 
0400  UT  on  December  7  (Day  341),  1982,  a  period  of  twelve  hours  Which  spans 
0800  to  2000,  Pacific  Standard  Time.  The  six  oblique- incidence  sounders  were 
active  during  this  time,  recording  ionograms  every  fifteen  minutes.  A  digital 
file  of  each  oblique  ionogram  was  processed,  by  Which  signals  in  excess  of  a 
threshold  value  D-0 . 15  were  plotted  in  oblique  ionogram  form,  and  sets  of  all 
these  plots  for  the  SSL-BCT  have  already  been  distributed  to  participants. 
These  plots  were  processed  to  obtain  true  height  profiles  for  Paths  (or  Zones) 
1,2,3,  and  3  during  1600-2200  UT  on  December  6.  It  was  inconvenient  to  use 
these  plots  for  Paths  4  and  6  since  they  frequently  did  not  show  an  E  layer, 
which  the  calculation  procedure  uses  to  convert  relative  time  delays  to 
absolute  time  delays.  A  peak-processing  method  was  developed  by  F.  Rhoads  of 
MRL  to  overcome  this  difficulty.  An  example  of  the  problem  and  this  solution 
to  it  are  shorn  in  Figure  22.  The  top  part  shows  the  oblique  ionogram  for 
Path  4  at  1800  which  resulted  from  the  D»0 . 15  threshold.  The  E  layer  is 
invisible,  and  there  is  not  much  high-ray  F  layer  in  the  trace.  The  result  of 
the  peak-processing  method  on  the  digital  file  is  shown  in  the  bottom  part  of 
Figure  22.  Some  E  layer  is  visible,  typically  somewhat  more  than  seen  here, 
and  the  high-ray  F  layer  trace  is  much  more  visible.  Ionograms  of  this  type 
were  processed  during  1600-1800  UT  for  Paths  4  and  6.  The  E  layer  parts  of 
the  calculated  true  height  profiles  are  not  expected  to  be  accurate  during 
this  time  for  Paths  4  and  6,  particularly  since  they  show  up  weakly,  but  the 
important  F  layer  parts  are  expected  to  be  reasonably  accurate  (cf.  Sections 
2C  and  2D). 

After  1800,  the  ionograms  processed  for  Paths  4  and  6,  provided  by  M. 
Daehler  of  NRL,  were  of  the  thresholded  type,  but  the  threshold  value  was  much 
lower.  The  1800  UT  ionogram  for  Path  4,  thresholded  at  D*0.02,  is  shown  in 
the  top  part  of  Figure  23.  The  E  layer  trace  is  strong  now,  for  our  purposes 
(cf.  Section  2D),  and  the  high-ray  F  layer  is  evident.  The  ionogram  processed 
at  a  later  time  is  shown  in  the  bottom  part  of  Figure  23.  Mote  the  appearance 
of  sporadic  E.  The  regular  E  layer  is  essentially  non-existent  at  this  time, 
since  the  solar  influence  is  practically  nil.  Note  also  the  nose  extension 


near  the  F  MUF.  It  Is  important  to  ignore  the  nose  extension  in  the  true 
height  profile  calculation. 

The  appearance  of  the  ionograms  for  the  paths  other  than  Path  4  is  now 
discussed.  For  the  most  part,  the  D-0 . 15  ionograms  for  Paths  1,2,3  and  5  are 
adequate  for  processing,  although  the  lower  threshold  ionograms  were  processed 
after  2200  UT  for  improved  accuracy.  Path  1  examples  are  shown  in  Figure  24. 
The  top  part  shows  a  D»0 . 15  ionogram,  and  D=0 . 06  for  the  bottom  part.  The  E 
MUF  region  is  just  visible  here,  but  often  it  isn't. 


Table  2  —  Location  of  midpath  and  vertical  ionosonde  points  in  Figure  2 


IDEMTIFICATIOH 

LAT.  (°m 

LOM.  (°M) 

Midpath  XI 

36.174 

119.67 

Midpath  X2 

36.498 

118.38 

Midpath  X3 

35.183 

116.97 

Midpath  X4 

38.673 

113.59 

Midpath  X5 

34 . 701 

119.44 

Midpath  X6 

41.896 

122.13 

Ft.  Ord  (R) 

36.652 

121.74 

LeMoore  (MAS) 

36.25 

119.95 

Pt.  Arguello 

Vandenberg  AFB  (VBG) 

34.600 

120.600 

Desert  Range  Station 
(DRS) 


38.667 


113.75 


A  similar  pair  of  examples  for  Path  2  is  shown  in  Figure  25.  Both 
examples  seem  to  show  a  sporadic  B  contribution  at  high  fraquancias ,  in 
addition  to  tha  regular  B  layar.  Figura  25  illustratas  an  apparent  problem 
with  tha  Path  2  data.  Thera  seams  to  ha  an  upward  equipment  drift  in  time 
delay,  superimposed  upon  tha  natural  one.  A  straight-edge  laid  along  tha 
sporadic  B  part  will  show  an  upward  tilt,  instead  of  tha  expected  flatness. 
This  drift  will  causa  some  inaccuracy  in  tha  true  height  profile  calculation. 
For  example,  the  calculated  values  of  hmB  and  haF2  are  expected  to  be  too  high 
from  this  effect,  but  the  results  are,  nevertheless,  expected  to  be 
qualitatively  accurate. 

Some  Path  3  ionogram  examples  are  shown  in  Figure  26.  The  tipper  part  is  a 
D»0 . 15  ionogram,  and  D-0.02  in  the  lower  part.  Vote  that  the  B  layer  is  faint 
and  incomplete,  and  not  very  much  of  the  high-ray  F  layer  is  seen  in  the  upper 
ionogram.  This  compromises  the  accuracy  of  the  true  height  profile 
calculations  for  Path  3  before  2200  UT,  although  the  trends  of  the  results  for 
the  Whole  group  of  ionograms  should  be  reliable.  After  2200  UT  the  results 
are  obtained  from  ionograms  like  the  bottom  one  in  Figure  26,  and  relatively 
good  accuracy  is  expected. 

Path  5  ionograms  are  similar  to  Path  2  ionograms,  as  expected,  since  the 
sounder  range  is  close  to  600  km  for  both.  Some  examples  are  shown  in  Figure 
27.  The  top  ionogram  was  measured  not  long  before  sunset.  The  bottom 
ionogram  was  measured  about  3.5  hours  later;  it  shows  pronounced  sporadic  B 
conditions.  Generally,  these  conditions  apply  toward  nighttime.  Relatively 
good  accuracy  is  expected  from  Path  5  ionograms. 

Ionograms  from  Path  6  are  shown  in  Figure  28.  The  peak-processing  method 
gives  the  upper  ionogram,  and  a  low  threshold  gives  the  lower  ionogram  a  short 
time  later.  Both  methods  result  in  ionograms  Which  should  yield  relatively 
good  accuracy  for  Path  6,  although  the  B  layer  is  somewhat  better  defined  in 
the  lower  ionogram.  The  question  arises,  in  connection  with  ionograms  like 
this,  about  which  part  of  the  thick,  low-ray  F  layer  trace  it  is  proper  to 
digitize.  It  is  presently  not  clear  Whether  to  digitize  the  leading  edge, 
middle,  or  trailing  edge  of  the  trace.  As  a  practical  matter,  however,  it 
doesn't  make  much  difference  in  our  scheme.  Where  the  observed  ionogram  time 
delays  are  relative  values.  We  choose  to  digitize  the  middle  portion;  the 
potential  loss  of  accuracy  is  relatively  insignificant.  As  stated  previously, 
the  magnetically  split  portion  of  the  high-ray  F  layer  portion  of  the  trace  is 
approximately  bisected  in  the  digitization  for  the  no-field  trace. 

3B.  True  Height  Profiles 

The  calculated  true  height  profile  results  for  Paths  1-6  in  Figure  2 
in  the  time  period  1600,  December  6  (Day  340)  to  0400,  December  7  (Day  341) 
are  given  in  Tables  3-8.  Bach  table  gives  a  time,  an  ionogram  number,  and 
either  nine  parameters  when  (fOB2,  hmB)  is  chosen  to  correspond  to  the 
maximum  plasma  frequency  calculated  in  the  B  layer  (cf.  Figure  4),  or  eleven 
parameters  When  (fOB2,  hmB)  is  chosen  to  correspond  to  the  first  reentrant 
point  (cf.  Figure  8).  In  the  eleven-perameter  case,  which  is  found  to  be 
prevalent  toward  nighttime,  the  additional  values  of  the  parameters  hsB  and 
fsB2  are  listed  directly  underneath  the  values  of  hmB  and  fOB2. 

_ Solar  Zenith  Angle  Dependence  and  Trend  Lines 

If  the  ionospheric  parameters  of  the  preceding  section  are  plotted 
vs.  the  time,  a  noisy  variation  about  a  background  diurnal  dependence  is 
obtained.  This  background  dependence  is  wave- like.  For  example,  the 
background  dependence  for  f0F22  is  seen  to  grow  during  the  morning  hours 
until  about  local  noon,  after  Which  the  dependence  falls  off  into  the 
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Fig.  26  —  Thresholded  oblique  ionograms:  D—0.15  (top)  and  D—0.02  (bottom) 
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Table  3  —  True  height  profile  calculations  —  Zone  1 


1555/1 
/i/a/9 
AZ9/9- 
/iio/f 

A6S/1 


1755/9 


t  M.6$  S.y///  /7S.95 
/  Mf.oj  4>./097 

3  /  0-3.26  6.3695  /7X'tO 

V  /69-17  1.165S  /7/-4/ 

r.  Jn£  4!  2.9713 /M  7 6 

1. /it,  / 1 

//.oil  /5°‘B*. 
/0.327  mu 
/0.2/2  /7S.76 
y/.O/i  /6o./o 
9X77/  176.5/ 
//■  /13  ns  77 
J*X*7  /61.31 
/ 0,097  /77.0a 
/*-96o  /72.30 
/0.7//0  /10.36 
1. 7(51/61.75 
729/9  779X9 
J/,550  /?%/* 
<7.  WU  /63. 
70.563  763-36 
7-VM  /97-7S 
9-  5/1*  /  99.71 
10,051  j99.il 
<1  A  (51 712,51 
IPm  1721  6 


hbE  =  fb 6*.  jth 

■hA  F  hntPl  •£< 

7X191  /9o.ir  3Z.5G  *19.73 
3.7J/55  2p7.  *7  4(2.166  *6j?6 
/0.917  2/3.11  32,133  25S./J 
/0, *60  220.1/  77.12/  273,73 
/0.VZ9  201,9/  7j.6tf  *66X9 
3.(2776  235.35  96.796  27253 
/0,?3f  / 13,21  H0.H71 2,69.79 
/0.9*2  113, St  96-739 262.23 
//.o3C  2/5.51  10, tr)  279-99 
7 0.163  171.  /!  *6.009  295.61 
/ 5335  20/Xl  97.093  *73.60 
2/.?96>  2 91,73  729.6!  2S5.12 
/?.  3 05  2*157  / ot.ot  27/21 
29-/66  /13.*6  60,27  *67.17 
/9.190  235.97  f  5,005  *12-00  |  ///J)S 
9/073  232.97  91.0/6  *19.92. 

/1: 937  2*9. 56  77.962  215.6/ 

7  7. 7/7  7/rjt  69. 69/  7.11.25 
22.9*7  207.7, Z  73.395  272.15 

H.  9563  ’*?.  17  2/.  531  &3.7* 

<7,397/  2?a,o9  19X5  *96.96 
/1.165  201V.  7/.  01 5  217.91 
/o.9/a  2*1.15  16.11*  217.90 
3.5176  2/S. 92  37129  *  *3. H 
9.*60(,  7X7.12  63.  925  73J3* 

I. 7255  23/., S3  77,117  299.7* 


Table  3  (Cont’d)  —  True  height  profile  calculations  —  Zone  1 
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Table  4  —  True  height  profile  calculations  —  Zone  2 
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True  height  profile  calculations  —  Zone  3 
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Table  5  (Cont’d)  —  True  height  profile  calculations  —  Zone  3 
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Table  6  —  True  height  profile  calculations  —  Zone  4 
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Table  7  —  True  height  profile  calculations  —  Zone  5 
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Table  7  (Cont’d)  —  True  height  profile  calculations  —  Zone  5 
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Table  8  (Cont’d)  —  True  height  profile  calculations  —  Zone  6 


night time  hours,  as  sxpsctad.  Further,  these  curves  ere  shifted  with  respect 
to  each  other  for  the  different  paths  in  a  manner  consistent  with  solar  zenith 
angle  variations.  The  sun  is  most  directly  overhead  at  different  universal 
times  (Greenwich  Meridian  Time)  for  the  six  midpath  points.  The  solar  zenith 
angle  X  at  an  observation  point  is  the  angle  made  by  the  line  of  sight  from 
the  observation  point  to  the  sun  with  the  vertical.  Which  is  directly  overhead 
from  the  observation  point.  From  the  observed  dependence  it  is  expected  that 
if  an  ionospheric  parameter  were  plotted  as  a  function  of  cos  x»  or  some  power 
of  it,  the  dependence  might  appear  as  a  noisy  variation  about  a  straight 
line.  This  straight  line  would  be  a  trend  line  from  which  a  prediction  of  the 
ionospheric  parameter  could  be  simply  made.  Evidently,  the  daytime  dependence 
on  solar  zenith  angle  is  more  fundamental  than  the  dependence  on  time. 

The  transformation  from  time  to  solar  zenith  angle  at  an  observation  point 
is  part  of  HIBIHUF  (Rose,  Martin,  and  Levine,  1978a],  a  computerized  algorithm 
for  predicting  F  MUF  values,  based  on  a  statistical  ionospheric  model  with 
sunspot  number  as  the  driving  parameter.  Later  updates  of  this  algorithm  are 
contained  in  [Rose  and  Martin,  1978b;  Rose,  1979;  Sailors  (1-6),  1984).  Part 
of  the  algorithm  has  been  used  to  convert  universal  time  to  cos  x  and 
cos  Xeff  at  F2  layer  heights.  Here,  X  is  the  actual  solar  zenith  angle 
and  X  eff  is  the  effective  solar  zenith  angle.  The  latter  incorporates  a 
correction  for  the  well-known  lag  between  the  solar  source  term  and  the 
ionospheric  response  in  the  upper  F  layer  region.  For  the  time  period  of 
interest,  the  lag  effect  is  small;  calculated  actual  and  effective  solar 
zenith  angles  are  nearly  equal.  The  plan  is  to  plot  ionospheric  parameters 
associated  with  the  upper  F  layer  region  (i.e.,  hmFl,  fOFl2,  hmF2,  f0f22, 
and  Y2)  as  a  function  of  cos  Xaff ,  or  some  power  of  it,  and  the  other 
parameters  as  a  function  of  cos  X,  or  some  power  of  it.  The  conversion  from 
time  to  actual  and  effective  solar  zenith  angles  is  given  for  midpath  points 
X1-X6  of  Figure  2  in  Tables  9-14. 

In  deciding  what  power  of  cos  X  to  use  as  the  independent  variable,  one 
may  recall  the  predictions  of  the  simple  Chapman  Layer  Theory  [e.g., 
Chamberlain,  1978],  whereby  the  rate  equation  for  electron  density  Mg  at  the 
maximum  of  the  layer.  Which  is  directly  proportional  to  fv  2  at  the 
maximum,  can  be  written  as 

dN 

-  C  +  q  cos  X  -  L  ,  (7) 


where  C  is  a  constant  source  term  which  is  not  dependent  on  x>  4  if  t  constant 
in  another  source  term  Which  gives  the  solar  influence,  and  L  is  a  loss  term 
which  has  one  of  two  common  forms: 


sN 


m 


0N_ 


( a  -layer) 
( 0 -layer) 


(8) 


In  quasi-equilibrium,  Which  applies  to  ionospheric  layers  most  of  the  time, 
the  left  side  of  (7)  can  be  approximately  set  to  zero.  The  result  is  that 
Mq  2  varies  linearly  with  cos  X  in  a -layer  theory,  and  Kg  varies 
linearly  with  cos  X  in  8 -layer  theory. 

Chapman-Layer  theory  is  too  simplified  for  the  upper  F  layer  region  of  the 
ionosphere,  where  other  source  terms  arise  from  diffusive  transport  along  and, 
to  a  lesser  extent,  across  magnetic  field  lines,  nevertheless,  it  will  be  of 
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0.516561 
0.5T3776 
0.509067 
0.502215 
0.4933ZT 
0.482438 
0,469594 
0. 454847 
0.438255 
0.419887 
o.3??65o 
0.354394 
0.329702 
~T>:  303673 
0.276413 
0.248032 

—ormsu 

0.188374 
0.15733 7 
“07175333“ 
0.093480 
0.060917 
“07078X08“ 
0,0123 02 
0.010856 
0,0<W>8& 
0.0074 62 
0.006585 
■  TrrooTiHi  i 
0.005128 
0,004525 


sms 


•o.v.v 
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Table  10  —  Time  to  solar  zenith  angle  —  Path  2  (cf.  Table  4) 


JHai. 

9 

XU  CooroLrt+.t<s  : 

L6T=  36.4979 

MO*  1?  DY«  6 

HOUR  (LIT) 

1.  OH  <  W  >  "  118.3753 

CQS(X) 

COS <  XFF  F ) 

16*004? 

pa/ 

0.292350 

0.281428 

%. 

16*  21)42 

390 

0.318583 

0.308141 

3 

16*5042 

0.34352? 

0.33360? 

V 

16.754? 

0.367064 

0 . 357706 

5" 

17.004? 

0.389114 

0 . 380357 

6 

17.2542 

0.409582 

0.401461 

7 

17.5042 

0.428385 

0.420933 

t 

17.7542 

0.445447 

0.438695 

<? 

18*004? 

0.460698 

0.454673 

to 

18*2542 

0.474077 

0.468803 

// 

18.7542 

0.495006 

0.491295 

13 

19.0042 

0.502472 

0,499567 

19.2500 

6.507821 

0. 505721 

•H 

19.5000 

0,511214 

0.509940 

•5 

19.7500 

0,512529 

0,312087 

/« 

20.0000 

0,511761 

oTslSISi 

n 

20.254? 

0.508847 

0,510083 

it 

20.5000 

0.503995 

0,506044 

/t 

20.7500 

0.497028 

0.499896 

no 

21.0000 

0.48804? 

0.491716 

U 

21.5000 

0.464162 

0.469401 

22 

21 . 7500 

0,449365 

0,455357 

23 

22.0000 

0.43274? 

0,439462 

a 

26 

22.2500 

0.414360 

0.421781 

22.5000 

23.0042 

6,394293 

0.349041 

?T.  462385 
0.35838? 

27 

23.254? 

0.324414 

0.334318 

29 

23.5042 

0,298468 

0.308893 

29 

23.7542 

p*y 

111 

0,271309 

0.28221 7 

30 

0.004? 

0.242760 

0.754075 

31 

0.2544 

0,213552 

0.225258 

32 

0.5042 

0, 183541 

0.195589 

33 

0.7542 

0.152749 

0. 165091 

3*/ 

1.0042 

0,121335 

O.Y33921 

35 

'  1.2544 

0.08939? 

0.102170 

36 

1.504? 

0.057156 

0.070074 

37 

1 «  754? 

♦  0I?4ASi 

0 . 037657 

39 

2.004? 

0.024651 

0.012646 

39 

2.254? 

0.0?4651 

0,011160 

Ho 

2.504? 

0.0?46M 

6,609849 

HI 

2.754? 

0.024651 

0,008691 

Hi 

3.004? 

0.024651 

0.007670 

HJ 

3.254? 

6.074651 

0.606/69 

hh 

3.504? 

0.024651 

0.005974 

Hs 

3.7542 

0,024651 

0.005272 

V* 

4.0042 

6.074651 

6.004652 

V7 

4.2542 

0.024651 

0,004106 
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Table  11  —  Time  to  solar  zenith  angle  —  Path  3  (cf.  Table  S) 
X3  Coore(tn^,ttS  • 

LAT~  35,1834  LON<W)«  116,9689 


LAT~  35,1834 
M8«  12  M- 
HOUR<UT) 


16*0853  0*  y 

16*3353  3  HO 

16*5853 


16*8353 

17,0853 

17,3353 _ 

17,5853 

17,8353 

18*0853 


18*3353 

18*8333 

19,0833 


1 , 3353 
1 , 5853 


3.0853 

3.3353 

3  7M5T 

3.8353 

4*0853 


c:ns<x> 


0 , 374808 
0,350867 
0.375516 


0.398654 

0.420190 

0.440036 

0.458113 

0.474348 

0.488676 


0,501039 

0.519622 

0,525841 


0, 104781 
0, 0731599 
0,040120 


0.006479 

0.006479 

6.606479 

0.006479 

0.006479 


COS(XFFK) 


0.313914 

0.340300 

0.365716 

0,389463 

0,411643 

0,432169 

0,450956 

0,467930 

0,483023 


0,496174 

0.516383 

0,523438 


19,3333 

0,529945 

—  ill  MINI 

19,5833 

0.531919 

0.531215 

19.8333 

0.531754 

0,531905 

20,0833 

0.529450 

20,3333 

0,525018 

0.526876 

20*5833 

0.518475 

0.521176 

20*8333 

—  II  II  1  II 

21.0833 

0.499169 

0.503521 

21,3333 

0.486484 

0.491637 

21 ‘,5833 - 

- 0,471843 - 

- 67477773 

21,8333 

0,455305 

0,461994 

22.0833 

0,436936 

0,444354 

23*0853 

0.346535 

0,356549 

23.3353 

0,320246 

0.330814 

23.5853 

0,292670 

0,303748 

23.8353  D. 

0,263916 

0,275460 

0,0853  Wi 

0,233905 

0,245839 

0,3353 

0,203224 

0,215530 

0.5853 

0.171725 

0, 184353 

0*8353 

0.139534 

0.152433 

0, 

0,086884 

0,053518 


A  4 

2.0853 

0,006479 

0,012170 

2.3353 

0,0 06479 

0.010740 

0,007382 

0,006514 

'077)65749 

0.005073 

0.004477 


Table  12  - 

Time  to  solar  zenith  angle  —  Path  4  (cf.  Table  6) 

X*4  CooroLn+fat  • 

ffo. 

1 

X 

^ 

LAT  =  .18* 
MCJ:=  1? 
HOUR(UT) 

673 

BY” 

L0N(M)»  113.5913 

6 

CHS ( X ) 

I'OStXFfFK  ) 

16,  0008 
16.2506 
16*5114 

pay 

3*0 

0.301289 

0.374639 

0.347594 

0,291490 

0,31534? 

0,338863 

V 

5 

6 

16.7506 
17.0006 

1 7 . 7506 

0,367268 

0.386343 

0.403819 

0,359090 

0,378777 

0,396896 

*r 

9 

1 

17.5006 

17.7506 

18.0000 

0,419623 

0,433689 

0,445972 

0.413371 

0.428134 

0.441138 

to 

II 

IX 

18.2508 

18.7500 

19.0000 

0,456393 

0.471513 

0.476173 

0,452300 

0.468946 

0.474388 

13 

iH 

IS 

19.7500 

19.5000 

19.7500 

0.478862 

0.479568 

0.478788 

0.477866 

0,479365 

0.478879 

1C 

;7 

it 

20 . 0000 
20,7500 
20.5000 

0,475027 

0,469800 

0,467677 

0.476410 

0.471969 

0,465573 

n 

XD 

XI 

20,7500 

21.0000 

21.2500 

0,453538 

0,442571 

0,479777 

0,457249 

0.447032 

0,434963 

22 

3-3 

xH 

21.5000 

21.7500 

22.0000 

0,415193 

0,398895 

0,380945 

0,421094 

0,405481 

0,388188 

xs 

26 

2-7 

22.2500 

22.5000 

23.0006 

0,361418 

0,340394 

0,294156 

0.369289 

0.348860 

0,303706 

2? 

2f 

50 

23.2506 

23.5006 

23.7508 

0,769185 

0,743099 

0,715976 

0,279218 

0,753574 

0 , 726850 

3/ 

32, 

33 

0.0008 

0.2506 

0.5006 

0*y 

31! 

0,187738 

0, 159079 
0,129728 

0, 198934 
0,170582 

0. 141490 

3H 

35 

36 

0.7506 

1,0006 

1 « 7506 

0,099840 

0,069537 

0,038946 

0,11181? 

0,081670 

0,051.189 

37 

39 

31 

~~  1.7506 
2,0006 
2,2508 

0.038946 

0,038946 

0,038946 

0,011231 

0,009911 

0,008745 

Ho 

Hi 

HI 

2,5008 

3.0008 

3.2506 

0,038946 

0.038946 

0.038946 

0,007718 

0,006011 

0,005305 

H3 

HH 

HS 

Hi 

3.5006 
3.7506 
4 . 0008 

0.038946 

0.038946 

0.038946 

0,00468? 

0,004132 

0.003646 

4.2506 

0,038946 

0.003218 

Table  13  —  Time  to  solar  zenith  angle  —  Path  5  (cf.  Table  7) 

X5  Coordinates  : 


LAT«=  34.70 07 
1?  OY* 
HOUR  (LIT) 


I.  ON  ( W )  •"  119.437 
COSOO 


17.5081 

17.8381 

18.0881 


18.3381 

18.8364 

19.0833 


19.3333 

19.5033 

19.8333 


20.0033 
20. 3333 


0,452926 

0.470527 

0,486244 


0,500014 

0.521441 

0,529003 


0,534553 

0.537962 

0.539217 


0.538312 

0.535251 


COS(XEFK) 


0.445232 

0.463562 

0.480035 


0,494506 

0,517632 

0.526020 


0.532419 

0.536685 

0.538802 


0,538760 

0,536562 


20,5833 

0.530047 

0,532214 

20,8333 

0,522720 

0,525734 

2 

l .  0833 

0,513299 

0,517150 

2 

1 , 3333 

0,501822 

0.506494 

2 

1 . 5833 

0.493809 

2 

1 . 8333 

0,472894 

0,479147 

22,0833 

0 .455558 

0,462566 

2: 

0,436398 

—  II  1  1  II  11 

22.5833 

0.415489 

0.423920 

22.8381 

0.392475 

0,401580 

23.0881 

0,368301 

0.378030 

23.3381 

0.342651 

0 , 352966 

23.5881 

0.315629 

0,326489 

23.8381 

0.287343 

0,298704 

i 

0.0881  P«.y 

0.257674 

0,269461 

1 

0.3381  59/ 

0,227280 

0.239475 

0.5881 

0. 195974 

0,208528 

I 

0.8381 

0, 163882 

0.176745 

1,0881 

0, 131132 

0.144753 

0,097857 

1.5881 

0,064189 

0.077667 

1.8381 

0.030264 

0,043839 

SV0W1 - 

0, 030,764 

“  0.013429 

2.3381 

0,030764 

0.011851 

2.5881 

0,030764 

0,010458 

2.8381 

0.030764 

0,009279 

3*0881 

0,030264 

0.008145 

3.3381 

0.030264 

0,007188 

3.5881 

0,030264 

0,006343 

3.8381 

0,030264 

0 « 005598 

4.0881 

0,030264 

0,004940 

Table  14  —  Time  to  solar  zenith  angle  —  Path  6  (cf.  Table  8) 

XC  Coe  £es  : 


4 1 . 8945 


MfJ=  12 
HOUR ( UT ) 


16.4278 

16*6700 


1*6700 
1*9200 
2 

2*4200 

2*6700 


U)N< W)~  t: 
COS(X) 


0 .252777 
0*274222 


1 248 


0.042279 
(>«  014378 


0.014378 
0* 014378 


COS(XFFF) 


92 

0.243542 

0.245400 


16.9200 

0.295204 

0.284861 

17.1700 

0.314928 

0,307092 

17.4200 

0.333310 

0,326014 

17.6700 

0,343546 

17.9200 

0.345741 

0.359415 

18.1700 

0.379451 

0.374152 

18.4200 

0.391944 

0,387094 

18.6667 

0.402436 

0.398247 

18.9167 

0.411346 

0.407864 

19.1467 

- 674TS531 - 

19.4167 

0.423937 

0.421852 

19.6467 

6.427524 

0,424162 

19.9167 

0.429290 

0.428657 

20.1667 

0.429226 

0 . 429324 

20.4167 

0. 427332 

0.428143 

20.6667 

0.423418 

- 6T435"174 

20.9167 

0.418699 

0.420378 

21.1667 

0.410798 

0.413789 

21.4167  0.461746  6" 465434 

21.6667 

0.390983 

0.395355 

21.9167 

0 . 378553 

0 . 383590 

22.1747 - 

- 0.344511 - 

- 67370T90 

22.4167 

0.348914 

0.355213 

22.9200 

0.313079 

0 . 320539 

23.1700 

0.293229 

0.301219 

23.4200 

0.272130 

0,280617 

23.6700 

0.249872 

0,258818 

23.-9760- ” 

- 0, 2245415 - 

072359T7 

0.1700 

0*4200 

0,201805 

0.211526 

0, 176730 

0.186791 

“077700 - 

- 0 , 156900 - 

0.161258 

0.9200 

0,124426 

0.135037 

1,1700 

0.097421 

0. 108240 

0.053375 
0.025541 

r 

0.009284 
0  <  008193 


3.4200 

0.014378 

0.005631 

3.6700 

0,014378 

0,004970 

4.1700 


0.014378 


0.003870 


'..n. 


some  interest  to  see  if  the  dependence  of  fOF2*  on  cos  X  eff  tends  to  be 
linear,  as  in  the  a  -layer  theory,  or  if  the  dependence  of  fOF22  on 
cos  x  «f f  is  aore  nearly  linear,  as  in  the  B -layer  theory.  In  addition  to 
these  choices,  the  MIHIMUF  ionospheric  model  asstnaes  that  fOF22  varies 
linearly  with  (cos  X  mff )^/2.  These  three  choices  have  been  investigated, 
as  shown  in  the  exaaple  of  Figure  29,  where  calculated  fOF22  values  for  Path 
1  in  the  tine  interval  1555-2210  (UT)  have  been  plotted  vs  (cos  X )1/2. 

These  values,  shown  as  pluses  in  Figure  29,  are  connected  by  straight  lines, 
and  the  result  is  a  noisy  variation  about  a  fairly  well-defined  background 
trend.  This  trend  is  generally  upward  with  increasing  values  of  the 
independent  variable.  It  is  believed  that  the  noise  present  arises  mainly 
from  wave-like  (e.g.,  TID's)  disturbances  in  the  ionosphere,  or  a 
superposition  thereof,  which  are  present  much  of  the  time.  Some  noisy 
appearance  may  result  from  calculational  inaccuracies,  which  are  difficult  to 
assess.  From  the  standpoint  of  making  predictions  of  fOF22,  it  is  the 
background  trend  in  Figure  29  which  is  of  interest  [Reilly,  1984b].  Three 
dotted  trend  lines  are  shown  in  Figure  29.  Trend  Line  1  is  a  linear 
regression  analysis  (best  linear  least-squares  fit)  of  the  calculated  f0f22 
values  on  (cos x eff )1/2,  which  corresponds  to  the  MIHIMUF  form.  Trend 
Line  2  shows  some  curvature  on  this  plot,  since  it  is  determined  from  the 
linear  regression  analysis  of  fOF22  on  cos  X  af f ,  which  corresponds  to 
the  B  -  layer  theory.  Trend  Line  3  also  shows  curvature,  since  it  is 
determined  from  the  linear  regression  analysis  of  fOF24  on  cosXeff,  Which 
corresponds  to  a -layer  theory.  Despite  the  different  appearance  of  the  three 
trend  lines  outside  the  region  of  calculations,  all  three  fit  the  calculations 
quite  well.  Still,  the  fit  of  Trend  Line  2  to  the  calculations  is  found  to  be 
somewhat  better  in  this  case,  and  in  the  fOF22  variations  for  the  five  other 
paths,  than  the  fit  of  Trend  Lines  1  and  3.  Indeed,  it  was  found  that  the 
extrapolations  based  on  Trend  Line  3  were  quite  erroneous  in  some  cases.  All 
the  other  ionospheric  parameters  were  treated  similarly.  The  general 
conclusion  was  that,  for  simplicity,  the  choice  of  independent  variable  as 
cos  X  or  cos  X ft  as  in  Trend  Line  2,  was  advisable,  since  good  fits  and 
reasonable  extrapolations  were  always  obtained.  For  the  display  of  results  of 
the  true  height  profile  calculations,  the  independent  variable  cosXs«f  will 
subsequently  be  chosen  for  the  parameters  hmFl,  fOFl2,  hmF2,  and  fOF22, 
and  the  corresponding  Trend  Line  calculations  will  involve  a  linear  regression 
analysis  of  these  parameters  on  this  independent  variable.  Similarly,  the 
independent  variable  will  be  cos  X  for  the  parameters  hmK,  fOE2,  hbF,  and 
fbF2.  The  difference  between  values  of  these  two  independent  variables  is 
very  small  for  the  time  period  of  interest,  but  it  can  be  significant  for 
summertime  operations. 

Of  further  note  in  Figure  29  is  the  predicted  behavior  of  fOF22  from 
MIHIMUF,  shown  for  the  observed  10.7  cm  radio  flux  value  of  210.4  for  December 
6,  1982.  This  translates  to  a  sunspot  number  of  167.  It  is  seen  that  the 
MIHIMUF  3.5  line  in  Figure  29  intersects  the  calculated  results  in  four 
places,  so  that  sounder-determined  fOF22  results  at  these  intersection 
points,  if  used  to  update  MIHIMUF  3.5  with  an  effective  sunspot  number 
(Uffelman,  Harnish,  and  Goodman,  1984],  would  determine  the  MIHIMUF  line  shown 
in  Figure  29.  This  line  is  seen  to  have  the  wrong  slope,  which  makes  its 
usefulness  for  predictions  questionable.  Although  MIHIMUF  is  apparently  based 
on  an  oversimplified  ionospheric  model  [Reilly,  1984b],  the  preceding  defect 
will  likely  be  shared  by  all  statistical  models  which  are  driven  by  only  a  few 
parameters,  such  as  the  sunspot  number.  It  will  be  argued  that  the  data  from 
a  configuration  of  sounders  can  be  analyzed  to  give  superior  capability  for  HF 
propagation  assessment. 
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Tables  3-14  can  b«  usad  to  plot  the  variation  of  tha  ionospheric 
parameters  with  aithar  cos  x  or  cos  X  ef f .  Thasa  plots  ara  displayad  by 
parameter  for  all  six  paths  in  tha  Appendices.  Appendix  B  displays  tha 
variation  of  fOF22  with  cos  x  ef f  for  X1-X6,  plus  tha  corresponding  results 
for  tha  vartical-incidanca  sounders  which  ware  available  for  comparison. 
Appendix  C  displays  tha  variation  of  hmF2  with  cos  Xaff  for  X1-X6 1  and 
Appandicas  D  and  B  show  this  variation  for  fOFl2  and  hmFl,  raspactivaly. 

The  variation  of  the  parameter  fbF2  with  cos  X  is  displayed  in  Appendix  F 
for  X1-X6 ,  and  corresponding  plots  for  hbF,  fOE2,  and  hmB  ara  given  in 
Appendices  G,  H,  and  I,  respectively.  The  results  are  discussed  next. 

The  results  for  fOF22  vs.  cos  X af f  are  plotted  in  Appendix  B.  They 
are  the  pluses  connected  by  tha  solid  lines.  Generally  speaking,  they  appear 
as  a  slightly  noisy  variation  about  a  roughly  linear  background  variation. 

The  results  start  out  in  the  middle  of  tha  figure  around  1600  UT  (0800  PST), 
advance  up  the  line  at  fifteen  minute  intervals  to  a  maximum  around  2000  UT 
and  come  back  down  the  line  until,  at  around  2200,  the  results  advance  up  to  a 
higher  line  and  conn  down  it  to  a  point  near  the  left  margin  of  the  figure, 
which  corresponds  to  the  beginning  of  nighttime  conditions  in  the  F2  layer  at 
around  0200  UT  (1800  PST).  From  here  (until  0400  UT)  the  values  tend  to 
decrease  along  a  steep  line  near  the  left  margin.  Calculated  results  are 
found  at  fifteen  minute  intervals.  They  are  seen  in  these  figures  to  be 
spread  out  more  for  larger  solar  zenith  angles.  Also  shown  in  each  figure  is 
a  set  of  dotted  trend  lines,  each  determined  by  a  linear  regression  analysis 
(best  linear  least-squares  fit)  on  a  portion  of  the  results.  For  example,  in 
the  results  for  XI,  four  such  trend  lines  are  shown.  Trend  line  1  (LI)  is 
found  by  fitting  points  1-9,  and  with  the  fifteen  minutes/point  conversion, 
this  corresponds  to  the  time  interval  1600-1800  UT.  Similarly,  L2  is 
associated  with  1600-2000,  L3  with  all  the  daytime  points  from  1600-0200,  and 
L4  with  2300-0200.  Greater  detail  of  the  correspondence  between  the  plotted 
point  and  time  can  be  found  from  Tables  3  and  9.  Close  inspection  reveals 
that  it  takes  1-2  hours  to  establish  a  trend  line  that  will  be  accurate  for 
predictions  ahead  for  one  hour  or  more.  Frequently  (except  for  X6),  the  trend 
line  established  between  1600  and  1800  UT  will  be  good  until  about  2200,  near 
which  there  will  be  a  transition  up  to  a  higher  trend  line.  This  trend  line 
can  be  established  between  2300  and  0030  UT,  and  it  will  be  good  until  about 
0200  in  this  ease.  For  best  accuracy  the  trend  line  should  be  updated  as 
often  as  possible,  but  it  is  convenient  that  predictive  capability  is  good  one 
or  more  hours  ahead.  This  allows  more  than  enough  time  for  the  data  to  be 
collected  and  processed  and  for  tha  trend  line  information  to  be  disseminated 
to  the  field.  Where  it  is  needed  for  the  SSL  or  other  HF  system  application. 

Available  vertical  ionogram  results  were  scaled  for  fOF22  values,  and 
these  results  are  interspersed  with  the  X1-X6  results  in  Appendix  B.  The  HAS 
results  were  placed  after  the  XI  results  for  comparison,  as  is  logical  from 
Figure  2.  Similarly,  VBG  was  put  after  X5.  There  was  only  one  DBS  ionogram, 
and  the  result  of  scaling  it  is  shown  in  the  X4  plot.  The  results  for  Fort 
Ord  (R  in  Figure  2)  are  shown  at  the  end  of  Appendix  B.  In  general,  the 
comparisons  between  calculated  f0F22  values  for  X1-X6  are  in  excellent 
agreement  with  corresponding  results  scaled  from  vertical  ionograms. 

Trend  line  information  for  the  eight  parameters  in  Appendices  B-I  is  given 
by  the  A  and  B  coefficients  of 
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Pj-AtPiH  B(Pi)  co*  X  (10) 

(P5-fbF2,  P6»hbF,  P7«fOf2,  Pg-hmE) 

Frequent  updating  of  the  A  and  B  coefficient*  ia  advisable  for  optinua 
accuracy  in  a  prediction  scheme,  but  not  necessarily  too  frequent  in  the  case 
for  f 0F22 ,  as  just  seen.  Other  parameters  may  require  different  treatment, 
as  trill  be  seen.  The  set  of  A  and  B  coefficients  for  the  ionospheric 
parameters  trill  be  of  use  to  SSL  and  other  HP  systems  in  the  field,  as 
discussed  later. 

The  effective  solar  zenith  angle  dependence  of  hmF2  for  points  X1-X6  in 
Figure  2  is  displayed  in  Appendix  C.  The  dependence  on  cos  X8ff  is  mildly 
increasing  in  the  morning  (local  time)  and  nearly  constant  in  the  afternoon. 

If  the  A  and  B  coefficients  are  frequently  updated,  there  should  be  no  problem 
in  using  the  associated  trend  lines  for  predictions.  If,  however,  the  trend 
lines  are  not  frequently  updated,  and  it  is  desired  to  use  them  for  long-range 
predictions,  special  procedures  are  indicated.  In  this  case,  the  results  of 
Appendix  C  suggest  that  if  a  trend  line  is  established  (using,  e.g.,  1-2  hours 
of  sounder  measurements)  during  early  morning  hours,  it  should  be  used  for 
predictions  only  during  the  morning  hours  (i.e..  When  cos  Xaff  is  increasing 
with  time) .  When  the  trend  line  is  established  during  the  early  afternoon 
hours  for  long-range  afternoon  predictions,  it  would  be  wise  to  assume  that 
B(hmP2)=>0,  in  Which  case  A(hmF2)  is  just  given  by  the  arithmetical  average  of 
the  hmF2  determinations.  Trend  lines  L5  in  the  plots  for  X4-X6  show  the 
pitfall  which  could  otherwise  develop. 

The  effective  solar  zenith  angle  dependence  of  fOFl2  is  shown  in 
Appendix  D.  Its  appearance  is  similar  to  that  for  fOF22.  The  parameter 
increases  roughly  linearly  with  cos  Xaff  in  the  morning  and  decreases,  again 
roughly,  linearly  with  this  parameter  in  the  afternoon  until  sunset,  beyond 
Which  point  the  decrease  is  steeper  in  appearance.  Frequent  updating  of  the 
trend  line  is  called  for,  but,  in  the  absence  of  this,  longer-range 
predictions  with  the  determined  trend  lines  have  to  be  made.  There  are 
limitations  to  the  success  of  this.  A  closer  look  at  the  dependence  of 
fOFl2  reveals  three  segments:  the  linear  increase  in  the  morning,  a 
transition  segment  in  the  early  afternoon  (e.g.,  between  2000  and  2200  UT), 
and  subsequently  a  different  linear  decrease  in  the  afternoon  until  sunset. 

The  trend  lines  established  in  any  one  of  these  segments  is  not  expected  to 
predict  well  in  the  other  segments.  This  applies  similarly  to  fOF22. 
Predictions  within  a  segment  should  be  reasonably  successful.  The  danger  of 
using  a  trend  line  in  the  transition  segment  to  predict  in  the  late  afternoon 
segment  is  shown  by  L5  in  the  X5  and  X6  plots. 

The  plot 8  for  hmFl  are  shown  in  Appendix  B.  The  background  dependence  of 
this  parameter  on  solar  zenith  angle  is  weak.  Frequent  updating  is  again 
advisable,  but,  in  lieu  of  this,  long-range  predictions  based  on  B(hmFl)*0 
should  be  quite  successful  and  safe.  If  this  is  not  assumed,  the  disparity 
between  LI  and  L2  for  X4,  and  between  L5  and  L4  for  X6  shows  what  could  happen 
in  long-range  predictions. 

The  solar  zenith  angle  dependence  of  fbF2  is  shown  in  Appendix  F.  The 
same  division  of  this  dependence  into  three  segments  (morning,  early 
afternoon,  late  afternoon)  is  seen  here  as  for  the  other  density  parameters, 
fOFl2  and  fOF22.  In  each  of  these  segments  the  background  dependence  of 
fbF2  is  increasing  with  cos  X  ,  but,  once  again,  the  trend  line  in  one 
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segment  cannot  bo  expected  to  necossarily  prod let  well  into  otbor  segments. 
Otherwise,  tha  trend  line  prediction*  should  work  quite  well. 

Plots  of  hbP  are  contained  in  Appendix  6.  This  height  parameter,  like  the 
others,  has  only  a  weak  dependence  on  solar  zenith  angle.  Frequent  updating 
of  the  trend  lines  will  give  optimum  accuracy,  but  safe  and  acceptable 
long-range  predictions  can  be  made  on  the  basis  of  the  assumption  that 
B(hbP)*0.  The  pitfall  of  making  long-range  predictions  without  this 
assumption  is  seen  from  L5  in  the  XS  plot. 

The  dependence  of  fOE2  on  solar  zenith  angle  is  displayed  in  Appendix 
H.  After  sunrise  this  parameter  exhibits  a  noisy  increase  up  a  line  until 
noon  When  it  starts  to  come  back  down  a  zigzag  path  about  a  similar  line.  In 
the  evening  it  begins  to  intersect  the  horizontal  axis,  and  in  some  cases  is 
spiky  in  appearance  at  this  point  as  sporadic  E  becomes  the  dominant 
contribution.  This  spikiness  has  to  do  with  the  treatment  of  the  K  layer  in 
the  true  height  profile  calculation.  Sometimes  it  is  handled  as  in  Figure  4, 
and  sometimes  as  in  Figure  8,  where  only  a  nominal  E  layer  contribution  is 
calculated.  The  X4  plot  shows  a  particularly  erratic  behavior  (e.g.,  LI)  in 
the  morning  region,  where  only  a  weak  portion  of  the  E  layer  was  visible  on 
the  ionograms.  The  spiky  appearance  toward  sunset  had  to  do  with  the 
ambiguity  of  Where  the  regular,  blanketing  E  layer  left  off  and  the  thin, 
patchy  sporadic  E  layer  began.  Frequent  updating  of  the  trend  line  is  called 
for.  For  long-range  predictions,  the  trend  line  determined  should  have 
positive  slope,  and  is  subject  to  the  three-segment  considerations  discussed 
previously,  although  prediction  from  the  first  segment  to  the  third  seems 
better  in  this  case  than  previously.  A  predicted  negative  value  of  fOE2 
should  be  set  to  zero. 

The  parameter  hmE  is  plotted  in  Appendix  I.  It  is  nearly  constant  until 
around  sunset  When  it  begins  to  jump  up  to  higher  values.  Again,  this  marks 
the  appearance  of  sporadic  E  layer  conditions,  and  the  conversion  from  Figure 
4  to  Figure  8  in  the  true  height  profile  calculations.  Frequent  updating  of 
the  trend  line  is  desired,  but  for  long-range  predictions  of  the  E  layer  and 
the  E-F  transition  region,  the  following  type  of  procedure  can  be  followed.  A 
trend  line  can  be  established  for  fOE2  in  the  afternoon  and  used  to  predict 
values  down  to  the  horizontal  axis.  Similarly,  a  trend  line  can  be 
established  for  hmE  during  the  same  afternoon  period.  It  should  have  the 
characteristic  that  B(hmE)«0,  or  else  this  condition  can  be  rigorously 
enforced,  to  be  safe.  Hence,  the  point  (fOE2,  hmE)  is  predicted,  and  the 
point  (fbF2,  hbP)  is  similarly  predicted  from  a  trend  line  determination. 

The  E-F  transition  line  from  the  former  point  to  the  latter  should  have 
positive  slope.  If  not,  move  (fOE2,  hmE)  up  the  reentrant  line  with 
negative  slope  until  the  new  E-F  transition  line  is  vertical.  This  will 
define  the  ionosphere  up  to  the  bottom  of  the  F  layer.  The  approximations 
involved  are  acceptable,  if  one  recalls  that  the  reentrant  region  is  ambiguous 
anyway,  the  E  layer  calculations  are  approximate,  and  the  important  F  layer 
calculations  are  insensitive  to  small  alterations  of  the  E  layer. 

It  has  been  estimated  in  the  preceding  discussion  of  the  plots  that  it 
might  take  one  to  two  hours  of  sounder  measurements  to  establish  a  trend  line 
for  the  background  dependence  of  an  ionospheric  parameter.  The  appearance  of 
the  parameter  is  noisy,  and  another  guideline  would  be  to  process  something 
like  a  full  cycle  of  noise  deviation  to  get  the  trend  line  information.  The 
recommendations  about  how  to  use  the  trend  lines  for  long-range  predictions 
were  based  on  analysis  of  data  during  a  twelve  hour  period  during  December 
6-7.  Vumerous  other  time  periods  remain  to  be  analyzed,  and  the  results  will 
be  of  interest  with  regard  to  the  preceding  and  further  observations.  For 
example,  the  nighttime  situation  has  not  been  fully  elucidated  by  the  present 
analysis. 
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3P.  Trend  Lines  Compared 

A  comparison  of  trend  lino  behavior  at  the  various  sounder  control 
points  trill  provide  some  insight  into  how  the  information  about  the  ionosphere 
at  sounder  control  points  (transmitter  positions  for  vertical-incidence 
sounders  and  midpath  points  for  oblique-incidence  sounders)  can  be  used  to 
infer  and  predict  the  ionosphere  along  unsounded  paths.  Recall  that  there 
were  basically  three  time  segments  in  the  daytime  during  each  of  Which 
ionospheric  trend  lines  could  be  established  and  used  successfully:  morning, 
early  afternoon,  and  late  afternoon.  Trend  lines  for  the  first  and  third 
segments  will  be  compared  here.  The  first  segment  is  given  by  the  trend  line 
L2  in  Appendices  B-I,  which  was  established  between  1600  and  2000  UT.  The 
third  segment  is  given  by  the  trend  line  L4,  Which  was  established  between 
about  2230  and  0200  UT,  except  for  fOB2  for  X2  and  X4,  where  early  afternoon 
hours  were  used.  The  trend  line  comparisons,  by  parasietar,  are  shown  in 
Appendix  J. 

The  L 2  plot  for  fOF22  in  Appendix  J  is  considered  first.  The  trend 
lines,  other  than  4  and  6,  are  close  to  each  other,  as  expected  by  the 
proximity  of  the  associated  control  points  in  Figure  2.  Further,  the  lines 
seem  to  occur  in  the  right  order;  8  and  0  are  close  to  1,  and  7  is  close  to 
5.  Curiously,  4  and  6  are  close  together  here.  One  might  try  to  justify  this 
on  the  basis  of  the  calculated  magnetic  dip  angles.  This  angle  is  measured 
upward  from  the  magnetic  field  vector  at  a  point  to  the  local  horizontal 
plane.  Values  for  X1-X6  are  found  to  be  61.7°,  62.2°,  61.3°,  65.0°, 

60.4°,  and  66.0°,  respectively.  The  formation  dynasties  of  the  upper  F 
layer  include  a  contribution  of  enhanced  diffusion  charge  transport  along 
magnetic  field  lines.  Hence,  in  the  vertical  distribution  of  ionospheric 
charge  and  the  parameters  of  the  F2  layer,  it  stands  to  reason  that  there 
could  be  a  dependence  on  magnetic  dip  angle.  The  L2  lines  do  seem  to  be 
strongly  correlated  to  the  magnetic  dip  angle  in  this  case.  In  the  next  plot, 
i.e,  the  L4  lines  of  fOF22,  it  is  found,  however,  that  this  correlation  with 
magnetic  dip  angle  is  not  as  impressive.  The  slopes  of  the  6  and  4  lines 
differ  substantially.  The  sounder  control  points  in  Figure  2,  Which  are 
geographically  in  close  proximity  to  each  other,  do  again  show  strongly 
correlated  trend  behavior. 

The  next  two  plots  show  the  trend  lines  for  hmF2.  The  L2  trend  lines  show 
that  this  parameter  tends  to  increase  in  the  morning,  and  the  L4  lines  show 
that  it  tends  to  stay  constant  in  the  afternoon.  Once  more,  the  different 
slopes  of  L2  lines  4  and  6  tend  to  indicate  that  geographical  proximity  of  an 
arbitrary  point  to  each  of  the  sounder  control  points  is  a  more  important 
consideration  for  a  spatial  interpolation  scheme  than  the  magnetic  dip  angle 
value  in  relation  to  these  values  for  the  sounder  control  points.  The  other 
trend  lines  of  these  plots  are  also  consistent  with  this  observation. 

The  trend  line  plots  for  fOFl2  are  qualitatively  similar  to  those  for 
f 0F22 .  The  same  consents  apply.  It  is  interesting  that  the  L2  lines  4  and 
6  are  fairly  close  to  each  other  again. 

The  trend  line  plots  for  hmFl  follow.  They  are  quite  similar  to  the  hmF2 
plots.  The  fact  that  the  L2  line  2  is  up  with  4  and  6  may  reflect  the  drift 
problem  mentioned  for  sounder  2  in  Section  2A,  which  would  tend  to  cause 
overestimated  values  of  the  height  parameters.  Similarly,  one  cannot  attach 
too  much  significance  to  the  decreasing  trend  of  L4.1ine  2  at  this  point. 

.  The  trend  lines  for  fbF2  Which  follow  tend  to  show  a  good  deal  more 
variation  than  seen  in  the  preceding  plots.  There  is  virtually  no  apparent 
dependence  on  magnetic  dip  angle,  nor  is  any  expected  this  low  in  the  F 
layer.  The  fact  is  that  the  lower  ionospheric  parameters  are  subject  more  to 
the  vicissitudes  of  the  true  height  profile  calculations,  including  the  B 


layer  approximations  and  a  naeassarily  approximate  reentrant  procedure.  This, 
and  the  fact  that  this  parameter  is  part  of  a  procedure  for  fitting  an 
idealized  profile  to  a  calculated  true  height  profile,  appear  to  account  for 
the  variability  in  the  fbF2  trend  lines,  nevertheless,  the  trend  line 
results  are  deemed  to  be  useful  in  the  interpolation  schema  to  be  discussed 
later.  There  really  is  not  very  much  variation  in  the  hbF  trend  lines  Which 
follow.  The  parameter  is  nearly  constant  on  a  given  trend  line.  The  last  two 
trend  line  plots  are  of  the  parameter  fOK2,  which  exhibit  the  strong  solar 
dependence  expected  during  the  daytime  period.  The  L4  lines  1  and  2  were 
determined  using  data  in  the  early  afternoon.  Hence,  they  do  not  include  data 
for  Which  sporadic  K  behavior  exerts  an  influence.  The  plots  for  hal  and  hbl 
have  not  been  included,  but,  disregarding  the  influence  of  sporadic  K,  they 
would  appear  to  have  only  a  very  weak  dependence  on  solar  zenith  angle.  In 
fact,  hbE»90  km  was  the  assumption  about  the  bottom  of  the  E  layer  under  these 
conditions . 

This  leads  into  the  question  of  the  specification  of  the  E  layer  under 
sporadic  E  conditions  and  nighttime  conditions.  Some  consideration  was  given 
to  this  in  Section  3C.  A  procedure  would  be  to  use  the  trend  lines 
established  in  the  afternoon  without  the  influence  of  sporadic  E  conditions 
for  hbE,  hmE,  and  fOE2.  Then,  as  previously  discussed,  the  fOE2  trend 
line  would  be  used  to  predict  this  parameter  down  to  the  horizontal  axis 
(f0E2*0).  The  hmE  trend  line  would  predict  hmE,  and  hbE*90  km  in  this 
case.  Then,  if  the  line  connecting  (fOE2,  hmE)  with  (fbF2,  hbF)  is 
reentrant,  i.e.,  has  negative  slope,  the  former  point  could  be  moved  up  a 
reentrant  line  (e.g.,  with  slope  -  /Z.5  km/MHz2)  toward  the  vertical  axis 
until  the  line  from  it  to  (fbF2,  hbF)  is  vertical.  Of  course,  the  other 
parameters  are  straightforwardly  predicted  from  their  trend  lines.  This  would 
complete  the  ionospheric  specification,  from  which  relatively  accurate  SSL 
calculations  could  be  made,  in  conjunction  with  the  interpolation  scheme 
discussed  next. 

4.  INTERPOLATION  OF  SOUNDER  RESULTS  TO  UNSOUNDED  PATHS 

In  SSL  a  measurement  of  frequency,  elevation  angle,  and  azimuth  angle  is 
made  on  the  incoming  wavefront.  These  are  initial  conditions  in  the 
calculation  of  the  raypath  to  the  HF  transmitter  through  an  ionosphere  Which 
is  specified  at  only  a  relatively  few  sounder  control  points.  The  necessary 
transition  from  ionospheric  knowledge  at  these  points  to  specification  of  the 
ionosphere  along  the  SSL  raypath  is  the  subject  of  this  section. 

As  is  well  known,  the  ionosphere  is  a  fluctuating  medium.  The  results  of 
Section  3C  are  further  evidence  of  this.  Rapid  fluctuations  will  often  not  be 
correlated  between  different  sounder  control  points  in  Figure  2,  but  the 
background  dependence  will  be.  Indeed,  statistical  models  of  the  ionosphere, 
though  often  inadequate  for  accurate  SSL  calculations,  make  global  predictions 
of  ionospheric  characteristics.  Which  represent  the  background  dependence. 
Ionospheric  fluctuations  affect  SSL  determinations.  Which  amount  to  fixing  the 
HF  transmitter  at  the  center  of  a  scatter  plot  of  fixes  (Reilly  and  Coran, 
1983b].  Fluctuations  affect  the  scatter  of  points,  but  it  is  hoped  that  the 
center  of  the  scatter  plot  is  associated  with  the  background  dependence  of  the 
ionosphere.  This  is  Where  the  trend  lines  of  the  last  section  come  in.  They 
are  assumed  to  represent  correlated  motion  of  the  ionosphere  on  the  time  scale 
used  to  establish  these  trend  lines,  and  are  found  to  predict  quite  well. 
Fluctuations  generally  tend  to  diminish  the  accuracy  of  the  SSL  determination, 
but  proper  use  of  these  trend  lines  should  represent  a  substantial  improvement 
in  SSL  technology.  This  remains  to  be  tested  in  future  research. 


The  specification  of  tha  lonoaphara  along  an  unsoundad  SSL  path  can  ba 
carrlad  out  by  intarpolating  tba  A  and  B  coefficient*  of  Equations  (9)  and 
(10)  in  Section  3C  for  tba  trend  lines  of  tba  ionospheric  parameters  at  tba 
sounder  control  points  (cf.  Figure  2).  An  example  is  shown  in  Figure  30, 
where  tha  sounder  control  points  of  Figure  2  are  shown  along  with  a  pair  of 
dashed  raypaths.  The  receiver  R  is  at  point  0  in  Figure  30,  and  the  control 
points  X1-X6  are  at  points  1-6  in  Figure  30.  Five  sample  raypath  points  (not 
on  straight  lines  in  the  general  case)  are  shown  as  Z-marks,  and  are  labelled 
a-e.  Only  the  trend  line  coefficients  for  1-6  have  been  found,  hut  it  will  be 
supposed  that  they  have  also  been  found  for  0.  Raypath  point  a  is  interior  to 
the  triangles  associated  with  the  points  0,1,  and  6  or  0,2,  and  6,  and  either 
set  of  three  control  points  is  thus  suitable  for  specification  of  the 
ionosphere  at  point  a  by  interpolation  of  the  trend  line  coefficients  at  these 
three  points.  Point  a  is  also  interior  to  triangle  0,4,6,  but  for  accuracy's 
sake,  only  the  smallest  possible  triangle  is  of  interest.  If  latitude  is 
denoted  by  L,  longitude  by  and  is  a  trend  line  coefficient  (A(P^)  or 
B(P£)  in  Equations  (9)  or  (10)),  then  a  suitable  interpolation  for  a  point 
interior  to  a  triangle  of  sounder  control  points  is 

Hi  -  AL  ♦  b£  *  C  ,  (11) 

where  L  and  j  refer  to  the  point  in  question,  and  the  constant  coefficients  A, 
B,  and  C  are  found  from  the  three  equations  Which  evaluate  equation  (11)  at 
the  three  sounder  control  points.  A  possible  alternative  scheme  is  to  perform 
a  weighted  average  of  values  at  the  sounder  control  points.  Where  the 
weighting  factor  for  each  sounder  control  point  term  is  the  inverse  of  the 
range  from  the  point  in  question  to  the  sounder  control  point.  A  calculation 
of  solar  zenith  angle  (effective  or  actual)  at  point  a  will  thus  enable  the 
specification  of  ionospheric  parameters  at  point  a  from  the  interpolated  trend 
line  coefficients  and  Equations  (9)  and  (10).  It  is  straightforward  to  also 
find  the  latitude  and  longitude  derivatives  of  these  parameters  at  point  a, 
from  Which  ionospheric  tilt  parameters  can  be  specified.  With  this 
information  the  next  raypath  increment  from  point  a  can  be  calculated  in  a 
ray- tracing  routine. 

Suitable  interpolation  triangles  for  raypath  points  b-d  are  246,  015,  and 
125,  respectively.  Point  e  is  just  outside  triangle  235,  but  it  is  close 
enough  that  extrapolation  from  these  points  with  (11)  should  be  pretty  good. 
Extrapolation  is  usually  more  risky,  however,  and  the  necessity  of  it  equates 
to  a  deficiency  in  the  sounder  network  distribution.  Hot*  that,  in  the 
absence  of  trend  line  coefficients  from  point  0,  the  triangle  for  point  a  is 
256,  and  point  c  would  involve  a  risky  extrapolation  from  triangle  125. 
Information  for  control  point  0  can  be  obtained  from  true  height  profile 
calculations  on  the  vertical- incidence  sounder  data  for  Fort  Ord.  The 
inclusion  of  the  control  point  at  the  receiver  is  generally  advantageous  for 
triangular  interpolation  possibilities. 

There  is  another  scheme  for  specifying  the  ionosphere  along  unsounded 
paths  which  may  be  advantageous  for  raypath  points  outside  interpolation 
triangles .  The  idea  would  be  to  use  a  statistical  ionospheric  model  Which 
incorporates  principles  of  ionospheric  dynamics  and  is  consistent  with 
averages  of  past  observed  ionospheric  behavior  on  a  large  geographical  scale. 
The  ionospheric  specification  at  sounder  control  points  could  be  used  to 
update  parameters  of  the  model  ionosphere  in  the  program,  which  would  then  be 
used  to  predict  the  ionosphere  along  the  unsounded  path. 

It  would  not  be  too  difficult,  in  principle,  to  incorporate  interpolation 
and  extrapolation  schemes  in  a  computer  ray-tracing  routine.  The  result  would 
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Fig.  30  —  SSL-BCT  sounder  control  points  (dots  -  cf.  Figure  2)  and 
sample  raypaths  (dashed)  with  sample  points  (X-marks)  along  them 
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be  «n  automatic,  relatively  accurate  SSL  capability.  Actual  Implementation  of 
these  Ideas  await  future  research  and  development. 

5.  MODIFICATION  OF  SSL  TBCHVIQUIS 
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A  generalisation  of  the  simple  SSL  method  of  the  Introduction  (Section  1), 
which  should  Improve  accuracy  without  sacrificing  too  much  in  the  way  of 
simplicity,  is  to  try  to  use  ionospheric  characteristics  at  the  SSL  midpath, 
if  it  is  believed  that  radiation  reached  the  receiver  by  a  one-hop  mode.  The 
first  step  is  to  use  ionospheric  characteristics  at  some  sounder  control 
point,  e.g.,  one  near  the  direction  of  the  incoming  wave  front,  in  the  manner 
discussed  in  the  Introduction.  It  was  found  that  the  vertical  ionogram  could 
be  simply  used  within  the  approximations  that  the  ionosphere  exhibits  no 
tilts,  i.e.,  it  is  spherically  symmetric,  and  that  the  ionosphere  is  nearly 
flat  for  the  portion  of  the  raypath  in  it.  A  not  too  much  more  complicated 
alternative  is  to  use  the  parameters  from  the  true  height  profile 
calculation.  As  noted  in  connection  with  Aquation  (4)  of  Section  2B,  the 
idealized  profile  fit  based  on  these  parameters  lends  itself  very  easily  to 
the  range  calculation  in  a  spherically  symmetric  ionosphere.  Having 
calculated  the  range,  it  is  a  simple  matter  to  infer  the  midpath  point,  and 
the  second  step  is  to  infer  the  ionosphere  at  this  point  and  repeat  the  range 
calculation.  Inferring  the  midpath  ionosphere  is  done  in  the  manner  described 
in  the  preceding  section.  It  involves  the  interpolation  of  ionospheric  trend 
line  parameters  at  other  sounder  control  points.  The  second  step  is  repeated, 
now  with  possible  inclusion  of  a  simple  tilt  model,  in  an  iterative  fashion 
until  the  calculated  HF  transmitter  position  converges  to  the  desired 
accuracy.  It  is  anticipated  that  only  a  few  iterations  would  be  needed. 

An  alternative  SSL  technique  is  more  complicated,  but  also  is  more 
accurate,  since  it  would  accurately  take  account  of  longitudinal  and 
transverse  ionospheric  tilts.  It  was  alluded  to  in  the  preceding  paragraph. 

It  would  involve  a  full-blown  ionospheric  ray-tracing  program,  with  initial 
conditions  given  by  the  SSL  measurements  and  with  ionospheric  characteristics 
updated  frequently  in  the  raypath  calculations  by  means  of  interpolation  of 
trend  line  parameters  for  ionospheric  parameters  determined  at  the  control 
points  of  a  sounder  network. 

Computer  automation  trill  facilitate  the  calculations  in  any  modification 
of  S&  techniques.  The  most  difficult  part  of  the  envisioned  SSL  system  is 
the  deployment  and  maintenance  of  a  sounder  network. 

6.  DISCUSSION 

In  the  Introduction  there  was  a  discussion  of  present  state-of-the-art  SSL 
technology,  in  Which  severe,  simplifying  approximations  are  made  for  the 
properties  of  the  ionosphere  between  the  receiver  and  the  HF  transmitter. 

Whose  location  is  to  be  determined.  The  oblique- incidence  and 
vertical- incidence  sounder  data  of  the  SSL-BCT  project  in  December,  1982, 
provided,  and  still  provide,  an  opportunity  to  assess  the  impact  of  proper 
accounting  of  ionospheric  properties  on  SSL  technology.  But  several 
computational  difficulties  stood  in  the  way.  The  analysis  hinged  on  the 
ability  to  process  the  ionograms  from  the  oblique- incidence  sounder  network, 
Which  constituted  the  major  part  of  the  data  base.  The  technology  had  been 
developed  for  vertical- incidence  sounders.  The  author  developed  a  large  part 
of  it  recently  for  oblique- incidence  sounders.  The  method  of  obtaining 
ionospheric  true  height  profiles  from  oblique  ionograms  is  delineated  in 
Section  2,  where  it  is  also  shown  that  a  true  height  profile  could  be 
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characterised  accurately  and  conveniently  by  a  sat  of  parameters.  The 
correspondence  between  the  calculated  true  height  profile  and  the  ionosphere 
at  the  sounder  midpath  is  also  confirmed  here.  The  method  was  applied  to  the 
set  of  oblique  ionograms  from  SSL-BCT  for  December  6-7  (UT) ,  which  spanned  the 
local  time  period  in  fifteen  minute  intervals  from  the  morning  hours  of 
December  6  after  sunrise  into  the  nighttime,  a  few  hours  after  sunset. 

Details  of  the  data  set  and  the  results  of  analysis  on  it  are  contained  in 
Section  3. 

It  was  found,  not  surprisingly,  that  the  daytime  ionosphere  was 
well-characterised  by  the  smith  angle  of  the  sun,  which  measures  at  a  point 
in  the  P  layer  the  angle  that  the  sun  makes  with  the  overhead  direction,  with 
some  accounting  for  the  time  lag  between  solar  perturbation  and  upper 
ionospheric  response.  When  the  ionospheric  parameters  were  plotted  against 
the  cosine  of  this  angle,  they  were  found  in  Section  3  to  show'  a  noisy 
behavior  about  quite  well-defined  trend  lines,  which  are  thus  used  to 
characterise  the  background  dependence  on  the  ionosphere.  These  trend  lines, 
updated  as  frequently  as  is  convenient,  could  be  used  to  predict  the 
background  dependence  of  the  ionosphere  at  the  sounder  control  points,  and  it 
was  shown  how  to  do  this.  It  was  also  seen  from  comparison  of  the  trend  lines 
at  the  sounder  control  points  that  the  background  dependence  at  these  points 
seemed  to  be  correlated. 

Correlation  of  the  background  dependence  of  the  ionosphere  at  sounder 
control  points  is  the  basis  for  specifying  the  ionosphere  on  raypaths  between 
the  SSL  receiver  and  the  HF  transmitter.  These  raypaths  practically  never 
include  the  sounder  control  points.  In  Section  4  it  was  shown  how  to  perform 
spatial  interpolation  of  trend  line  parameters  at  the  vertices  of  sounder 
control  point  triangles  for  the  purpose  of  specifying  the  ionosphere  along 
raypath  portions  within  these  triangles.  Other  possibilities  for  inferring 
ionospheric  properties  along  unsounded  paths  were  mentioned.  These  involve 
the  updating  of  the  parameters  of  a  statistical  model  ionosphere  with  sounder 
results,  and  this  method  may  prove  to  be  especially  useful  for  inferring 
ionospheric  properties  on  raypaths  Which  are  outside  the  aforementioned 
control  point  triangles. 

In  Section  5  it  was  discussed  in  detail  how  a  proper  HF  propagation 
assessment  with  ionospheric  sounders  will  result  in  a  modification  of  SSL 
technology.  The  data  from  either  vertical-incidence  or  oblique-incidence 
sounders  can  be  used  for  this  purpose  to  an  essentially  equivalent  effect. 

The  important  thing  is  to  establish  a  strategic  distribution  of  sounder 
control  points.  Practical  considerations,  such  as  cost,  safety,  and  security 
of  transmitter  and  receiver  location,  will  also  be  used  to  determine  the  mix 
of  oblique-  and  vertical- incidence  sounders.  Trend  line  parameters  at  the 
sounder  control  points  can  be  radioed  to  the  SSL  receiver  position.  Where  the 
HP  transmitter  location  can  be  computed.  An  improvement  in  accuracy  over 
present  methods  can  be  obtained  by  inferring  the  ionosphere  at  midpath  between 
the  receiver  and  transmitter,  in  the  manner  discussed  in  Section  5.  Of 
course,  if  the  raypath  between  receiver  and  transmitter  is  inferred  to  involve 
more  than  one  hop,  the  ionosphere  should  be  specified  at  all  of  the 
ionospheric  reflection  points.  With  this  information,  standard  approximations 
can  be  used  to  calculate  the  HP  transmitter  location.  For  greatest  accuracy, 
a  ray-tracing  routine  with  ionospheric  specification  along  the  raypath  can  be 
used  in  the  manner  discussed,  lo  formidable  obstacles  appear  to  stand  in  the 
way  of  these  modifications  to  SSL  technology,  once  the  commitment  to  deploy 
the  sounder  network  has  been  made. 

Several  tasks  remain  in  the  analysis  of  SSL-BCT  data.  Several  of  the 
other  time  periods  in  the  data  base  can  be  analyzed  to  provide  additional 


information.  In  particular,  techniques  for  characterizing  and  interpolating 
nighttime  ionospheric  properties  need  to  be  developed.  Further,  the  actual 
computation  of  HF  transmitter  positions  can  be  carried  out  using  the 
methodology  of  this  report.  The  appropriate  incorporation  of  this  into  a 
ray-tracing  program  needs  to  be  further  developed.  Greater  accuracy  in  the 
interpolation  of  ionospheric  parameters  trill  be  obtained  from  additionally 
processing  the  vertical  ionograms  at  Fort  Ord.  The  success  of  the  various  SSL 
techniques  can  be  quantitatively  evaluated  by  comparing  actual  and  calculated 
HF  transmitter  positions.  Basic  information  about  ionospheric  effects  will  be 
obtained.  Comparison  of  these  results  with  the  contractor  results  will  also 
be  of  interest  in  this  regard.  Further,  the  ionospheric  effects  have  to  be 
subtracted  out  to  get  a  better  understanding  of  errors  in  measurement 
techniques  used  by  the  contractors. 
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APPENDIX  A 


Comparison  Between  Calculated  and  Experimental  Vertical  Ionoarama 

Path  4  oblique  ionograms  (shown)  are  analyzed  to  obtain  true  height 
profiles  (nine-parameter  solutions),  which  are  used  to  calculate  equivalent 
vertical  ionograms.  These  are  compared  in  the  following  with  ordinary  mode 
traces  of  vertical  ionograms  measured  at  Port  Ord  (FO)  and  Desert  Research 
Station  (DRS).  See  the  description  of  the  text.  Times  for  the  associated 
ionograms  are  listed  in  the  legends  for  these  comparisons . 


I 


r 

i 

j 


n 

\ 


r 

r 


SSL  0IS#2  ZONE  4  ERIE, CO  TO  FT  ORD.CR 


1982  348  16:00:02  ZULU 


4  6  3  10  12  14  16  18  20  22  24  26  28  30 


FREQUENCY  < MHz ) 

OIS#2  ZONE  4  ERIE, CO  TO  FT  ORD.CR 


1982  348  16: 15:02  ZULU 


10  12  14  16  18  20  22  24  26  28  30 
FREQUENCY  (MHz) 


SSL  0IS#2  ZONE  4  ERIE. CO  TO  FT  ORD.CR 


2  4  6  8  10  12  14  16  28  20  22  24  26  28  30 


FREQUENCY  (MHz) 


2  4  6  8  10  12  14  16  18  20  22  24  26  28  30 


FREQUENCY  (MHz) 


V 


K  N  ^ 


H  S  0 

oU 


+  *  ■ 

+•  *  a 

4-  x  ■ 


Si 

♦  * 

♦  • 

:  < 

- i"T 

♦♦  "I 

♦  *1 

x» 

* 


g 

&  P 

i>  iu 


o 

<9 

CD 

CD 

CD 

CD  CD 

<9 

<9 

CD 

CD 

in 

CD 

in 

CD 

in  CD 

in 

CD 

in 

M> 

ur> 

in 

rt 

it 

ro  ro 

CM 

CM 

•H 

SkCE 

®  CD 


VERTICAL  IONOGRAM 


o  o  ^ 

O  O  ‘O 
Oo  00 


K,  J  o 
O  §  ^ 


i  *  ■ 

X  » 

+  X  a 


.  I 
%  ♦ 


© 

CD 

© 

•3» 

CD 

© 

in 

© 

in 

© 

in 

© 

in 

-t 

ro 

ro 

w 

<N) 

OP  ® 

®  in 


84 


VERTICAL  IONOGRAM 


TIME  DELAY  mSEC  ~SIG  dBm  i  TIME  DELAY  mSEC  ~SIG  dBm 


» 


ZHM 


UERTICAL  10N0GRAN  +  OX  oX'/S.'OX 

Xx  x  PRS  oZ.15 


TIME  DELRY  mSEC  -SIG  dBm  •  TIME  DELAY  mSEC  -SIG  dBm 


SSL  0IS#2  ZONE  4  ERIE, CO  TO  FT  ORD.Cfi 
1982  352  15:30:02  ZULU 


2  4  6  8  10  12  14  16  18  20  22  24  26  28  30 

FREQUENCY  (MHz) 


SSL  OIS#2  ZONE  4  ERIE, CO  TO  FT  ORD.CR 
1982  352  15:45:02  ZULU 


TREND  LII 


P  3*/0j/6oo  -  ,0*/Oo  (uT ) 

4U  LINES* 

PTS.1,9 
P7S. 1.16 


XtH 


*6 ;  D3H< 

TREND  Lll 


Ft.  Orel  CR 

TRENO  LII 
l;  ;  PTS 
L2:  PTS 


APPENDIX  C 


Plots  of  haP2  vs .  cos  x  aff 


115 


XS :  D3H0  ?  v6oo  -  D  3H  t  >  oho o  cut') 

TREND  LINES: 

LI*  PTS.1,9 
L2*  PTS.1.16 


''COS' 


I 

APPENDIX  B 


eff 


Plots  of  hnFl  vs.  eosx 


i3 io  > /&oo  -  P3  )  oqoo  Cor) 

LINES: 


<-OS  Xe£f 


*5:  Ddyo>Ke>oo~  03H t >  0H00  C  ut} 

TREND  LINES: 


X6  :  0  3HlO>/6OO  -  D3H!  )  OHOO  iUT') 

TREND  LINES: 

LI:  PTS.1,9 
L2:  PTS .1,17 


''CO  s 


CO  s 


'CO  S 


",■*  ■>  ■;*  -r* 1  :  *  v»  *  •  '.“■  ^  \r- i”*  ."•  ‘ ■*.  .■%"  v v"  ^  w  v. *im .»  a-  *—■ 


\ 

.s. 


I 


►V 


APPENDIX  6 

Plot«  of  hbF  va.  cog  x 


i 


S 


£ 


i  ' 

►y 


t':: 

•  • 


£ 


143 


.%  /•. 


Cj-0)oo/rO  '  th£Q  -  009/  (0h£<]  \1X 


9*0 


/COS 


SO-> 


—  P  3  V  /  >  o  Cut') 

10  LINES > 

PTS.1,9 


/COS 


APPENDIX  I 


Plot 3  of  hmg  vs.  cos  X 
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Trend  lines  for  each  of  seven  ionospheric  parameters  are  plotted  together 
for  all  the  sounder  control  points.  Curves  1-6  refer  to  control  points  X1-X6 
in  Figure  2.  Additional  trend  lines  for  fOF22  from  vertical-incidence 
sounders  are  plotted  for  comparison.  These  are  identified  in  the  accompanying 
legend  for  the  plot.  For  each  parameter  there  are  two  plots,  one  for  the 
morning  trend  lines,  identified  as  L2  in  Appendices  B-H,  and  one  for  the  late 
afternoon  trend  lines,  identified  as  L4  in  Appendices  B-H.  The  lines  are 
drawn  in  the  intervals  for  which  they  were  determined  by  linear  regression 
analysis . 
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APPENDIX  K 

Comparison  of  Calculated  and  Exact  Profiles  in  the  SSL-BCT 
Prototype  Calculations 


From  a  nine-parameter  idealized  profile  with  a  parabolic  E  layer,  one-hop 
oblique  ionograms  are  calculated  for  the  ranges  of  the  SSL-BCT  sounder  network 
(cf.  Figure  2).  The  true  height  profile  calculation  procedure  is  carried  out 
on  these  ionograms,  and  the  nine-parameter  solutions  are  plotted  as  X-marks  in 
the  ensuing  figures.  The  exact  model  profile  is  shown  by  the  solid  line 
curves,  for  comparison.  Section  2C  gives  a  full  discussion  of  the  procedure 
and  results. 
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